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Our Society 


We wish to draw members’ attention to three changes introduced this month 
on the page giving particulars of the Society and its Officers, changes which we are 
sure will be warmly welcomed by all. 

Firstly, the offices of President and Vice-President replace those of Chairman 
and Vice-Chairman of Council. Secondly, the Scottish and Yorkshire Provisional 
Branches have been granted full Branch status. Both changes are referred to more 
fully in the News Section, on p. 387. 

Finally, the following statements have been introduced to amplify the title of the 
Society: 

Society for the Promotion of the Science and Technology of Aerospace and Astronautics 


A Founder Member-Society of the International Astronautical Federation 


The first subtitle indicates that the Society is concerned with both scientific and 
engineering aspects of our field, and with both the atmospheric fringe (aerospace) 
and outer space itself. It should not be taken as implying that we are uninterested 
in the legal, economic and sociological problems involved—perusal of past issues 
of the Journal will show that this is not the case. The second statement indicates 
our participation in the work of the appropriate international scientific union— 
the I.A.F. The Society played a leading role in the formation of the I.A.F., drew 
up its original draft Constitution and has since been a major contributor to its 
work and finance. 

No change in the Society’s policy is implied by the first subtitle. Neither is it the 
prelude to a change in name. From time to time we receive suggestions for 
changes, including the dropping of the word British, or the changing of the word 
Interplanetary to Rocket, Aerospace, Astronautics, Space, etc. However, the title 
of the Society was chosen with care and will be of value for many years yet—it 
would remain valid even were interstellar travel to become a reality. It is un- 
necessary for us to keep changing our name to accord with the current catchphrase 
of the moment, and our interests are not restricted to a particular form of propulsion 
—indeed, propulsion is but one of the many subjects we study. 

The need for retaining the word British seems obvious, for ours is a Society 
founded, incorporated and having its headquarters in the United Kingdom and 
although we have a large international membership we do not wish to encroach, 
or even appear to encroach, upon the territories of other astronautical societies. 
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A word on our relation to other British societies is not out of place here. We 
are not disturbed to their interesting themselves in appropriate aspects of space- 
flight. On the contrary, we welcome and encourage such interest, and the series 
of joint symposia which we have arranged has this as one of its objects. By 
arousing interest elsewhere we are promoting the science and technology of astro- 
nautics. Even where the Society does not directly participate in meetings of other 
bodies, we often assist them by providing or suggesting speakers. 

There is no need to draw up lines of demarcation, for the activities of various 
societies must inevitably overlap. Thus, although this Society is not concerned 
with aircraft as such, we cannot avoid dealing with winged re-entry vehicles. 
Similarly, although it would be inappropriate for the Royal Aeronautical Society 
to spend the majority of its time discussing interplanetary travel, it quite rightly 
covers such topics as guided missiles, the use of rockets for assisted take-off, and 
aerospace activities. 

Covering as we do all possible extra-terrestrial actions of man and his instruments 
(plus terrestrial construction, launching and support facilities) it is inevitable that 
our interests will have some common ground with those of many other organiza- 
tions. However, ours is the only society in Britian dealing in the broadest sense 
with all aspects of aerospace and astronautics. 

In disccusing the astronautical interests of other professional societies it is 
appropriate to mention the Royal Society. This eminent body has a unique place 
in the history of science, being the oldest and most distinguished of all learned 
societies. Now, in its tercentenary year, by its appreciation of the vast possibilities 
of space research, it has demonstrated that it is still in the forefront of scientific 
thought. The B.I.S. hopes that the United Kingdom Government will display 
the same degree of enlightenment in taking the necessary steps to ensure that the 
nation shall play an active role in exploring this great new frontier. 


Recovery from Orbit 


We offer our congratulations to the members of the American and Russian teams 
whose work on the Discovery and ‘Arknik’ satellites respectively has been crowned 
with success. These first recoveries from orbit are great technical achievements. 

As we go to press little information has been made available about the Russian 
spaceship. However, it seems evident that it will not be long before man himself 
will orbit the globe and re-enter the atmosphere. Indeed we would not be sur- 
prised if one of the Russian delegates to the 1961 International Astronautical 
Congress were able to report his personal observations of the other side of the 
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GRAPHICAL PLOTTING OF INTERPLANETARY ORBITS* 


By M. VERTREGT,} Fellow 


ABSTRACT 


A method is described in which a conic section can be represented by straight lines ; it enables a point on a conic section 
to be constructed easily if the length of the radius vector, or the angular pericentre distance is given. 


Methods are given for the construction of tangents to 


conic sections, for the graphical construction of the time needed 


for a body to move from one point to another along an elliptical orbit, and for the plotting of an interplanetary orbit and of 


the direction and magnitude of the transfer impulses. 


I. THE BASIS OF THE CONSTRUCTION 


FROM an astronautical standpoint, most graphical con- 
structions of conic sections have the disadvantage that 
they start from the centre, whereas what we need is a 
focal construction. Moreover, it is mostly not possible 
to construct a special point, e.g., a point with a given 
radius-vector or with a given pericentre-distance. The 
usual procedure is to construct several points, to connect 
these with a smooth line and thus to find the desired 
point indirectly. 
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The focal equation for a conic in polar co-ordinates 
(Fig. 1) is: 


r = pi(l + €cos¢) (1) 
We have: 
cos ¢ = — x/r 
from which: 
r=piex .. i Pa (2) 


This equation suggests a simple construction. 

Let us assume that the parameter p and the eccen- 
tricity « of a conic section are given (Fig. 2). Then, 
from a given point F, the focal point, we erect a per- 
pendicular FH with length p. 





* Manuscript received 29 May, 1959. 
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We draw a line MHD in such a manner that tan « = e. 
Then every perpendicular RS with distance x from F 
will have a length r in accordance with Equation (2). 

Now, if we.*race a circle with centre F and radius r, 
which intersects RS in Q, then Q is a point of the ellipse. 
To complete the diagram we draw two lines FC and FD 
which make angles of 45° with the base MA. From 
points C and D we drop perpendiculars CP and DA on 
the axis MA. 

Now if 
CP = PF = — x, 
Equation (2) gives 
x, = — pl + @). 
But this is the pericentre-radius of the ellipse. There- 
fore point P is the pericentre. 

Similarly, DA = x, = pil — e), which is the 
apocentre-radius, and therefore A is the apocentre of 
the ellipse. 

If we construct the points G and £, symmetrically to 
the major axis PA (Fig. 3), we get a trapezium, which 
we can consider as a transformation of the enclosed 
ellipse, because it shows many characteristics of the 
ellipse. These are: 

(1) The line FH is the parameter p of the ellipse. 

(2) The distance FP is the pericentre-radius of the 

ellipse. 

(3) The distance FA is the apocentre-radius of the 

ellipse. 

(4) PA is major axis of the ellipse. 

(5) OX is the semi-major axis of the ellipse (O is the 
centre). 


+ Mozartlaan 73, The Hague, Netherlands. 
351 


1959-60, Vol. 17 








352 Vertregt: Graphical Plotting of Interplanetary Orbits 




















D 
a 
x 
V lo 
Cc 
Pp 
a 
M P F OS A 
G 
E 
Fic. 3. 


(6) The slope of the line CD, or tan «, is equal to the 
eccentricity € of the ellipse. 

(7) The perpendicular MV in M is the directrix of the 
ellipse. For if MV is the directrix and VQ || MS, 
then according to a well-known definition of a 
conic section: FQ = RS = MS tana = « VQ. 


Il. CONSTRUCTIONS FOR ELLIPSES, 
PARABOLAE AND HYPERBOLAE 


Elements of a conic section are: 
(1) The semi-major axis a. 

(2) The semi-minor axis b. 

(3) The parameter p. 

(4) The eccentricity «. 


These elements are connected by the relations: 
b=av(l—e), 
p = a(l—e*). 
If two of these four elements are known, the trapezium 
CDEG can be constructed. 








1. ELLIPse 


If the length of a radius r of an ellipse is given (AK, 
Fig. 4), we draw RK||PA, draw RS | PA, and FQ = RS, 
then Q is a point of the ellipse and ¢ is the pericentre- 
angle of the radius-vector. 

Alternately, if the angle ¢ is given we proceed as follows 
(Fig. 5). Draw a line FL, which makes an angle ¢ with 
PA, set off AN = FL. Draw FN, which intersects CD in 
R. Draw RS | PA, which intersects FL in Q. Then 
Q is a point of the ellipse and the length of the radius- 
vector is FO = RS. 


2. PARABOLA 


For a parabola (Fig. 6), we have to know only the 
parameter p. Therefore we erect a perpendicular FH 
with the given length p on MS. We circle MF = FH 
and draw the line MR. Bisect MF and draw a per- 
pendicular PC in P. Then P is the pericentre of the 
parabola. As the eccentricity of a parabola is 1, then 
a = 45°. 


3. HYPERBOLA 


The construction for a hyperbola is as follows (Fig. 7): 
If the axis a and b of the hyperbola are given, we set off 
a length a and a perpendicular to this of length 6. The 
hypotenuse of the triangle is: 

(a? j 5?)! e .. [a? a a (e? —_ 1)}} S on de. 


Now from point F we set off FO = ae and from O, OP = 
a,andOZ =a. Oisthecentre ofthe hyperbola. Draw 
a perpendicular PC = PF and draw a straight line 
through Z and C, which intersects the horizontal in M. 
For a hyperbola, «>1, thus «>45°. Make PW = b, 
than OW is the asymptote of the hyperbola. For a 
parabola or a hyperbola the figure is not a proper 
trapezium, because it is open at one side. 
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CONSTRUCTION OF THE TANGENT 


To draw a tangent in a point Q of a conic section we 
first derive the equation of the tangent (Fig. 8). 
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then: 
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(PF — ye 
[(p + ex? — ep? 
[e(p + ex) — x] 





dx ~ [(p+ex?— x}? 


(re — x) 


~ (F— x2 


= tany 





(4) 


We now construct the tangent as follows: 

Draw a line through P and Q which intersects the per- 
pendicular on O in T. Draw TU\|PA. Then WQU is 
the tangent in Q. 


This can be demonstrated as follows: 

TO = QS.PO/PS, and WS = SA.QS/(TO — QS), 
tan y = QS/WS = (TO — QS)/SA 
((QS.PO/PS) — QS)/SA 
= QS (PO — PS)\(PS.SA). 


Now: 


OS = y = (r° — x1, 
PO — PS = FO — FS = ae — x = [pe/(1 — &)] — x 
= (re — x)(1 — €'), 
PS = [pP(1 + 0))+x=["F -—ex)(1+9)+~x 
=(r+ x)(1 + 2), 
SA = [pl — 2] — x = [(r—ex)(1— 8] — x 
= (r — x)/((1 — €). 


Thus: 





até —x/[r+xr—-x 
tan y = (# — ap. 3 /[* 7 a 
= (re — x)/(r? — x*)'”. 


in accordance with Equation (4). 
For point H, x = O, and tan y = «. 
CD is a tangent to the ellipse in point H. 


Thus the line 
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Fig. 9 gives the construction of the tangent at point @ 
of an ellipse on the other side of the centre. A similar 
construction can be used for a hyperbola. A hyperbola 
has two branches and the centre O (Fig. 7a) is situated 
in the middle of the line connecting the two foci. Call 
the pericentre of the first branch P,, and that of the 
second branch P,. Draw a line from Q to P,, which 
intersects the perpendicular through O in a point T. 
From T draw a line parallel to the major axis, which 
intersects a perpendicular through P, in a point U. 
Draw a line through Q and U; this line is the tangent 
in Q. 
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For a parabola, « = 1, and thus according to Equation 
(4): 
OS/(FW + SF) = tan y = (r — x)/(r? — yp? 
= (FF — x) P(r + x) = vir + x). 
As SF = x, FW must be equal tor. Therefore the con- 
struction of the tangent is as follows (Fig. 10): Make 
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FW = rand join W with Q, then W@Q is the tangent in Q. 
In this system a circle is represented by a square 
and its radius is, of course, half the side of the square. 
To draw a tangent at point Q (Fig. 11), prolong FQ 
till it intersects CG in U, then make WF = UF. WQ 
is the tangent in Q@. The demonstration follows from 
the equality of the triangles WQF and UPF. 
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IV. INTERPLANETARY ORBITS 


The intersection of an ellipse and a circle is represented 
in Fig. 12, and the intersection of an ellipse with two 
concentric circles in Fig. 13. This is the case of an inter- 
planetary orbit. 
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As in a previous paper, I assume that the orbits of the 
planets are circular and coplanar. If the orbits of two 
planets, say Earth and Venus, are represented by the 
squares C,D,E,G, and C,D,E,G,, then every straight 
line which intersects the lines C,D, and C,D, represents 
the s/ope-line of an interplanetary orbit. This orbit 
can be elliptical, parabolic or hyperbolic. To construct 
such an interplanetary orbit we prolong line R,R, till it 
intersects the line FC, in C to the left and FD, in D to the 
right of F. When the slope of the line R,R, is equal to 
or greater than 45°, this line will obviously not intersect 
FD, and we get a parabolic or a hyperbolic orbit. 

Construction of the two points of intersection Q, and 
Q, of the ellipse with the two circles gives us the angle ¢ 
or the angular distance of the voyage. y is the angular 
distance of the planets at the moment of departure of the 
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spaceship from Earth. (Q, is the situation of the planet 
of destination at the moment of departure; Q, is the 
situation of the Earth at this moment and Q, is the point 
of encounter of the planet and the spaceship. 

When the line R,R;, intersects the lines C,D, and C,D, 
in the points C, and D, we get a so-called Hohmann- 
orbit (Fig. 14). Inspection of this figure shows that: 
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a = (r, + 13)/2. 


If: 
q = a/r, and n = r,/r, 
then: 
q=(l+na)/2.. “ S% (5) 
and: 
e=(,—r)/(1 +r) =U—va(l+n) .. (6) 
and furthermore: 
FH = p = 2nr,/(1 + n) re “te (7) 


V. THE DURATION OF A VOYAGE 


The duration of the voyage following an elliptic orbit 
is represented by the equation!: 


T = a, ({M)2 ={(n —<esar) —(— esin7,)}, 
where a is the semi major axis of the ellipse, 


fis the gravitational constant, and 
M is the mass of the central body (the Sun). 


For a complete revolution of the Earth around the Sun, 
we get: 


T' (= 1 year) = r°’?. (f{M)“?? 27, 


where r is the radius of the Earth’s orbit, and, 
M is the mass of the Sun. 


Then 
T (in years) = a*, ({M)~, { (, — € sin r,)— 


— (r%, — esin 7a) } [P,({M)"2, 2n 
=(¢° */2m).{ (r — € sin 7,)— 
— (7, — € sin r)} (8) 
where : 
q is the semi major axis of the elliptic orbit in A.U. 
7, = arc cos [(1 — r,/a)/e] 2% A (9) 
T, = arc cos [(1 — r,/a)/e] a .. i 


where r, and r, are radii of the ellipse and a is the semi 
major axis. 


If we call: 

™ — €sint, = T; ia bie we GN 
and: 

T, —€SiNt, = T; Rs - -» (12) 
then Equation (8) becomes : 

T = (q°/2/2)(T; — T2) be o 


The calculation of this form is rather cumbersome and a 
graphical construction will save much time, if the 
required accuracy is not too high. 
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To achieve this, we construct a common cycloid 
(Fig. 15) of which the equations are: 


x = 4(27 — sin 2r) Te ao. 
y=+(1 —cos2r) .. i sp. Tee 
and: 
dy — 
Fe = tr wi ‘4 a ae 
The equation of a tangent to the cycloid is: 
pagel FF Te 
y y= a & x) - vtec Se 
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where y’ and x’ are the co-ordinates of the tangent and y 
and x the co-ordinates of the tangential point. 
Introducing Equations (14), (15) and (16) into Equa- 
tion (17) we get: 
x’ =[(y’ —D/cotr] +7... .. (18) 
We now make: 


y =rla ie .. (19) 
and construct : 
7 = arc cos [(l —r/a)/e] .. .. (20) 
which gives : 
ria = 1 — ecost 
Then: 


y =1—ecosr 
Substituting this in Equation (18) we get: 

x’ = — [ecost/cotr] +7=—=7r—esinr=T’ (21) 
The graphical construction of 7’ goes as follows: 

We mark a point r/a on the y-axis (point A, Fig. 15). 
From point A we describe a circle with radius «, which 
intersects the line: y = 1 (line BC) in D. We draw a 
line AE through A and D. We now draw a tangent to 
the cycloid FG||AE. Then draw AF parallel and FH 
perpendicular to the x-axis. Then: 

OH = x’ =7r—esint=T’ 
The construction is very easy if we use squared paper 
and parallel rules. 
If we write the equation: 
T, = @2.7;/22 

in logarithmic form, we get: 

log 7, = 3 log q — log 2m + log T; a> 
We can solve this equation by means of a special slide- 
rule (Fig. 16). First we transform the linear scale A of 


the cycloidal nomogram into a logarithmic scale 7”. 
The slide has two logarithmic scales for g, which have 
a relative length of 1-5 times scale T’. The fixed scale 
T, which corresponds with the undermost scale of the 
slide, is shifted to the right over a length of log 27 
w:th regard to scale T’. 


Example. 


Assume that a body is following an elliptical orbit around 
the Sun, the elements of the ellipse being a = 0-625 A.U., 
« = 0-645. It is required to construct the period of time, in 
years, needed for the body to proceed from a point of radius 
r, = 1-0 to the perihelion. 

We have r,/a = 1-0/0-625 = 1-6. 

Using the cycloid (Fig. 16), and knowing r,/a and «, we 
construct 7’ on scale A as described above, and find the value 
T’ = 2:52. We now set the slide in such a manner that the 
mark 1-0 of the uppermost scale of the slide coincides with the 
value 2-52 on scale 7’ (as shown). 

Since gq = a = 0-625, we now determine the point on the 
lowest scale T coinciding with a point 0-625 on scale g. This 
point is 0-199, which is the time in years it takes the body to 
reach perihelion. 


In this manner it is possible to determine the time T 
very rapidly, and, if the cycloid is drawn on a sufficiently 
large scale and the logarithmic scales are of sufficient 
length, an accuracy of about 0-3% is easily attained. 
The construction is less exact when the angle 7 is nearly 
0°, 90° or 180°. Note that for r/a = 1, T’ = m/2 —e«. 


VI. PLOTTING AN ORBIT 


A common task in astronautics will be to plot an orbit 
to a certain planet, starting from Earth at a certain 
moment. Let us assume that the radius of the planet’s 
orbit is n A.U., and that the angle which the Earth and 
the planet are making at the moment of departure is #% 
(Fig. 13). If T is the time in which the spaceship must 
reach the planet, then: 
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T = (q°/27).[(7, — € sin7,) — (rT. — sin 72)] 
= (¢°"/27).(T; — T)) 
Now the time which the spaceship, travelling from Q, to 
Q., takes to reach the planet, must, of course, be equal 
to the time which the planet takes to travel from Q, to 


Q>. 
Thus: 


T =(n®?/2n\d+y) ..  .. (23) 


where ¢ is the angular distance of the voyage. 
Then: 


¢=2nTn32?—y% .. .. (24) 


The value of the form 277.n-*/* can be found by the use 
of the logarithmic scales. 

To find an orbit which satisfies this condition we try 
several pairs of radii r, and r, for which the angular 
distance is ¢. To achieve this we assume r, to make an 
angle ¢, (Fig. 13) with the major axis, then: 

ti = $2 + ¢. 
We now construct the points R, and R, and trace the 


slope-line CD. It is not necessary to trace the other 
lines. From this line we measure g and e€: 


q =(DD' + CC’)/2 Dralion <a 
« = (DD' — CC’)\(DD' + CC’) .. (26) 


From this we derive the time 7, using the nomogram. 
We do this for three or four values of ¢, and set off the 
resulting 7’s against the ¢’s in a graph. 

We can determine which value of ¢, tallies with the 
required value of 7. Then for this value of ¢, we can 
construct the slope-line and find the values of q and e. 


VU. THE DIRECTION AND THE 
MAGNITUDE OF THE IMPULSE 


Now, as the elements of the ellipse are determined, we 
must find the direction and the magnitude of the impulse, 
which must be given to the spaceship for it to follow the 
prescribed orbit. These also can be constructed easily. 

If we call the elliptical velocity of the spaceship at the 
moment of departure V, and the circular velocity of the 
Earth V,, then: 


V3 = p(2/r, — 1/a) = V3 2 — n/a) 
where »p = fM. 


Thus: 
V,/V, = (2 — r,/a)* ane ae! 


Now, in Fig. 17, the focal point of the ellipse is F, the 








radius of the orbit at the moment of departure is FO = 
r,, and W@Q is the tangent at the ellipse in point Q. WF 
is the direction of the major axis of the ellipse. 

Then, if A is the length of the perpendicular from F 


on the tangent, 
h? = WF* — WZ? = [(y/tan y) — x}? — A*/tan® y. 

From which: 

h = (y/tan y— x)/((1 + I/tan*y)}* .. -- (28) 
Now according to Equation (4): 

tan y = (r,¢ — x)/y, and y = (r,2 — x*), 
and, derived from Equation (2): 
x = (r, — p)e. 

By substituting these relations in Equation (28) we find 
after some transformation : 


h? = r,p/[2 — r,(1 — *)/p] -- (29) 
And, as: 
Pre (1 €*), 
we find: 
h? = r,p/(2 — r,/a) =~ .. (30) 


Combining Equations (27) and (30), we find: 
V,|Ve = (Vp) ih 


Starting from this equation, V, can be easily constructed, 
Fig. 18. 











Vie 
Fic. 18. 


We prolong r, = FQ with a length p (FH in Fig. 13), 
bisect (r, + p) and draw a circle with radius (r, + p)/2 
from the centre M. In Q we erect a perpendicular on 
FQ; the length of this line up to the point of intersection 
with the circle is Vr,p. 

Now we choose a length on the perpendicular which 
represents the velocity V,. If we make V, = 29.8 cm., 
then all velocities have the scale of 1 cm. per km./sec. 
The direction of the velocity V, is the tangent QW to the 
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ellipse in point Q. By constructing a parallelogram of 
velocities we find the direction of the velocity V,, and 
the magnitude of the impulse, expressed as a velocity, 
which we must give the spaceship in order to adapt its 
velocity to the velocity on the elliptical orbit on the point 
of departure. 

At the point of arrival our spaceship must adapt its 
velocity to the velocity of the planet of destination. The 
circular velocity of this planet is: 

VS = blr, = w(n.r,) = V2in, 
or: 
V/V, = 1/Vn. 


This value can be constructed as shown in Fig. 19, which 
needs no elucidation. 











Then the velocity V,,, necessary to adapt the velocity 
of the spaceship to the velocity of the planet can be con- 
tructed in the same manner as the velocity V,,. 


LIST OF SYMBOLS 


semi major axis of a conic section. 
semi minor axis of a conic section. 
gravitational constant. 
perpendicular from the focus on a tangent. 
radius of the orbit of a planet (A.U.). 
parameter of a conic section. 
semi major axis of a conic section (A.U.). 
radius of a conic section. 
y rectangular coordinates. 
mass of central body (the Sun). 
duration of voyage along elliptic orbit. 
velocity of spaceship. 
circular velocity of Earth. 
circular velocity of planet of destination. 
a V2, transfer velocities. 
angle between tangent and major axis. 
eccentricity of a conic section. 


= fM. 
= arc cos[(1 — r/a)/e]. 
angular distance to pericentre of a point on a conic 
section. 


angular distance of a voyage. 
angular distance between Earth and planet at moment 
of departure. 


& & wet OS Sees TS SS OU SR 
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NOTE ON “ESTABLISHMENT OF A CIRCULAR SATELLITE ORBIT 
BY DOUBLE IMPULSE”* 


By Professor G. LEITMANN,? Ph.D., Fellow 


ABSTRACT 


In an earlier paper’ conditions were deduced for the existence of a relative minimum value of propellent expended in a 
two-impulse transfer into a circular satellite orbit. In this note the relation of relative and absolute minima is investigated. It 
is shown that in all cases the relative minimum is also the absolute one. 


I. INTRODUCTION 
AN eartier paper’ dealt with the determination of the 
direction of launch and the staging ratio leading to 
minimum propellent consumption for a_ two-stage, 
impulsive rocket required to carry a given payload into 
a circular orbit of prescribed radius. Launch from rest 
and from the rotating Earth were discussed. Conditions 


yielding a stationary (relative) minimum value of 
propellent consumption were determined and found to 
depend critically on the ratio of the specific impulses of 
the two stages. 

If a function possesses corners, i.e., discontinuous 
derivatives, and if the independent variables are bounded, 
the absolute extrema may occur at corners or at the 
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boundary of the domain. Hence, it is not clear a priori 
that the relative minimum found in reference (1) is the 
absolute minimum. It is the purpose of this note to 
investigate whether or not the absolute minimum value 
of fuel consumed is indeed the relative minimum value 
previously deduced. The nomenclature of this note is 
the same as that of reference (1). 


Il. COMPARISON OF MINIMUM VALUES 
The function to be minimized is: 


F=1vvit Vi+ v2— 
1 


iy ey ay 4+ © “stg 
: 7, VPA V2) 4 z ar 
where 
V2+(V,+ ve =p(A—Vv3) .. (2) 


For real values of velocity the arguments of the square 
roots in Equation (1) cannot pass through zero in the 
interior of the domain. Hence F possesses no corners 
and we need investigate only its values at the boundary 
of the domain. From Equation (2) it is seen that the 
independent variables chosen, V, and V,, are bounded 
by 

V2+ pV? = pA i ss (3) 
so that 
0<V2<pA, 0O< V2<4A ‘a (4) 


Since F possesses only one extremum with respect to 
V,, at V, = 0 —and it was shown to be a minimum— 
we need concern ourselves with the values of F only at 


V,.=OandV,=0,V,=+VA .. (5) 


Of these values of F the one at V, = V,=0 is an 


extremum so that only the one at V,=0, V, = VA 
(the negative value of V, is of no interest in physically 
realistic problems) must be found. It is 


intr =t 1 <a 
ga ya TA 4 Vi ig: Ty 
Cy Ce 


F possesses two extrema and for 


Qh n\? Ve 
—— =>} — eee ee 7 
C2 Te (3) p(V pA — vz) ” 


the minimum occurs at 
V,.=V,=0 “4 yd (8) 


and is given by 
— 1 = Cy " l — 
F= > [ ve ( rs “) r| + = Vreg, (9) 
In comparing the values of F given by Equations (6) 
and (9), respectively, we ask if 


Vi+Az vPA ( -25) 9, .. (0) 
2 


To answer this question we must choose values for 
C,/¢C, and v, which make the right-hand side of Inequality 
(10) largest. With respect to c,/c, this means choice of the 
equality sign in Inequality (7). In that case both sides 
of Inequality (10) are positive so that we may square 
them. Thus, we have in place of Inequality (10): 


2x 
> aa ll — , ee ee 
12 (1 — p*(l — x) pase: at (11) 
where 
nh en a ae 
i Ma 7, Vite 


For any p, the right-hand side of Inequality (11) takes 
on its largest value when 


x=p .. $a cin ae 


We may note that this corresponds to actual values of 
v, Only for sufficiently small values of p. With Equation 
(12), Inequality (11) becomes 


mt Bed 
121 i_-- rs . €e 
Thus, for p < 1 the right-hand side of Inequality (13) is 
less than its left-hand side, i.e., even for the “worst” 
choice of parameters the relative minimum is less than 
the value of F at the boundary, hence is absolute. For 
start from rest this result follows directly from Inequality 
(10), of course. 

If Inequality (7) holds reversed, the relative minimum 
precedes the boundary, hence cannot exceed the value 
at the boundary and is absolute. 
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A NOTE ON RELAY SATELLITES* 
By S. W. GREENWOOD, + B.Sc., M.Eng., A.M.I.Mech.E., A.F.R.Ae.S., Member of Council 


ABSTRACT 


Clarke’s proposal’ for the establishment of “stationary” relay satellites in the 24-hr. circular equatorial orbit is extended 
to the other planets of the solar system. General expressions for the orbit radius and the area of surface visible from the 


satellite are derived and values for each planet tabulated. 


INTEREST has long been expressed in the 24-hr. circular 
orbit around the Earth in its equatorial plane, partly 
because of the novelty of the prospect of a free body 
rotating around the Earth at the same angular velocity 
as that of the Earth about its own axis, and partly 
because of the practical possibilities of using such a 
satellite as a relay station for radio transmission over the 
portion of the Earth’s surface visible from it.! 

In contemplating the extension of this principle to 
other planets in the solar system, it is interesting to con- 
sider whether any new aspects of the problem arise. 
The dynamics of the artificial satellite’s motion would 
clearly be the same—it would be necessary to move in 
the equatorial plane of the planet at the same angular 
velocity as that of the planet about its axis, and hence 
in a circular path about the planet. A significant 
difference in viewpoint, however, is that the setting up 
of the satellites will almost certainly have preceded 
landings on the planet concerned. It seems to the 


associated with the selection of the type of orbit under 
consideration : 

(a) Continuous observation of a selected area of the 

planet over a prolonged period would be possible. 
Such an area might be of interest either for scien- 
tific reasons, or because it was being considered 
as a possible area for landing and preliminary 
exploration. 
In the period following a landing, the satellite 
could maintain continuous contact with the 
exploring party visually and by radio, and could 
act as a relay station for contact between ground 
parties widely separated from one another. 


These suggestions are of a general character, and will 
obviously be dependent upon sufficiently favourable 
circumstances in any given case for their being put into 


practice. 
Data for the planets of the solar system are presented 


(d) 













































































author that the following specific advantages would be in Table I. Caution should be exercised in connection 
TABLE I 
| R _A_ Notes 
___Body eae: ack t, hr. mit Ro, miles | a R, miles R, Amax | (see below) 
Mercury 88 x 24 1550 0-36 | 197,000 127 0-99 
Venus 13 x 24 3850 | 087 | 136,000 35-3 097 CY 1 
Earth es oe ee 26,200 6-63 oss 
Mars. .| (24-62 210 «|| 2038 | 12,660 6-03 0-83 | 
Jupiter .| 986 | 44350 | 6 | 100,500 2-27 056 2 
Saturn | 10-23 2% 37,550 | 13 | 69,300 1-85 046 ~~ CO 2 
Uranus... ..|:1082,s«|—=S 4,650 | 107 37,800 | 258 0-61 2 
Neptune... 14. hei Eide de ee | 47,400 | 3-42 0-71 2 
Pluto. | 64 x 24 | 2000 | rh 45,600 | 228 096 | 3 





LIST OF SYMBOLS 


time for one rotation in orbit, hours. 


gravitational acceleration at the planet’s surface. 


area of surface visible from satellite. 


Ry a radius of planet, miles. 
a = Bol(¥o)wartn 

& = 

R = radius of orbit, miles. 
A om 

Amax oe 


area of one hemisphere of the planet’s surface. 


NOTES: 


“Giants” do not rotate as solid bodies. The rotational period is given for the equatorial region. 


1. Rotational period very uncertain. 


3. All data very uncertain. 


Surface gravitational acceleration estimated assuming planet has same mean density as the Earth. 
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with much of the material included in the table, owing 
to the unreliability of some planetary data. The 
figures are, therefore, approximate and will have to be 
revised and the derived data calculated more accurately 
at a later date. 
The radius of the orbit under consideration is given 
by: 
R32 = 44-7 tR, a’? +. (1) 


The fraction of the area of a hemisphere of the planet's 
surface that may be seen from the satellite is given by: 


A/Amax = 1—R,/R ee oe (2) 


These relationships are derived in the Appendix. 

The results of the calculations also are presented in 
Table I. The planets appear to form three distinct 
groups for the purposes of this study. 


Group 1: Mercury, Venus and Pluto. 


The relay satellite orbits are at large distances from 
these planets, mainly owing to their long rotational 
periods. The effect of this is to permit almost a complete 
hemisphere to be viewed in each case, but at such a large 
distance that the value of doing it is questionable. 


Group 2: Earth and Mars. 


The orbits of relay satellites for these planets are 
relatively closer, and a relatively large proportion, over 
80°, of a hemisphere of the planet is visible. The use 
of such satellites appears quite practicable here. Clarke’ 
has proposed three spaced at 120° to each other for sub- 


stantially complete radio coverage of the Earth. 


Group 3: Jupiter, Saturn, Uranus and Neptune. 


All the “Giants” rotate rapidly, and the associated 
relay satellite orbits lie fairly close to the planets. This 
makes the satellites excellent observation stations, but 
reduces the extent of the area visible to about half a 
hemisphere. The orbits lie within the orbits of the closest 
known moons of the planets, but it may be that small 
moons do in fact exist in the regions of interest which 
could be used as bases. In the case of Saturn, the orbit 
lies within the ring system and it is debatable whether 
such an arrangement would be practicable. 

If a leisurely and more extensive visual inspection of 
the planet is contemplated in place of a prolonged 
examination of a particular area, this might be accom- 
plished by moving in an orbit at a slightly different radius 
from the one considered here. By a suitable selection 
of orbit, the planet will appear to the observer to rotate 
at any desired rate within limits. 

Studies similar to the above were carried out for some 
of the larger moons in the solar system, but without 
success. The influence of the gravitational field of the 


planet concerned would make the establishment of a 
circular orbit around the moon at the distances involved 
impossible without the application of power. 
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APPENDIX 
1. Assumption 


The orbit of the satellite is only affected by the gravita- 
tional field of the planet around which it rotates and is 
circular. 


2. Derivation of Equation | 


Putting ¢ = time for one revolution in orbit 
R = radius of orbit 
R, = radius of planet 
2 9 = gravitational acceleration at the planet's 
surface 
Qn R32 
Then t= Rog? 
, t Ry g,'” 
(ger = cat 
Rearranging, R vn 


For purposes of computation, with R and R, meas ured 
in miles and ¢ in hr., this may be written 
R3/2 — 44-7 tR,a/? 
where a = go/(g,)kartn 


3. Derivation of Equation 2 
Area of surface visible from satellite, 
B 
A= | 2» Resin ada 
0 
Where § is the complement of the angle between lines 
drawn from the satellite (a) to the centre of the planet 


and (6) tangential to the surface of the planet respectively. 
The integral becomes 


A = 2a R32 (1—R,/R) 
The maximum possible area visible is the area of a 
hemisphere, 
Amax = 2a Ri 
Thus A/Amax = 1—R,/R 
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OPTIMUM FREQUENCY SELECTION FOR INTERPLANETARY COMMUNICATION* 


By Professor SURESH CHANDRA,?+ M.Sc. 


ABSTRACT 


The power required for one-way radio communication from Earth to various planetary bodies has been calculated 
at various frequencies by choosing suitable minimum detectable power, molecular atmospheric absorption, receiver band- 
width and antenna gains. It has been found that Jupiter is not easily accessible to radio communication, owing to its large 
distance and heavy absorption by ammonia present in its atmosphere. However, as discussed in the paper, short-interval 
contacts may be established at about 100 Mc./sec. under most favourable conditions. A frequency of about 3000 Mc./sec. 
has been recommended for communication with the Moon, Mars and Venus. Under the conditions assumed in this paper 
the transmitter power required for communication with the Moon, Mars, Venus and Jupiter at 3000 Mc./sec. using con- 
ventional receivers is 7 W., 230 kW., 70 kW. and 8-5 MW. respectively. The power requirement for the transmitter 
decreases if a receiver employing a maser preamplifier is used. The power of the transmitter at 3000 Mc./sec. calculated 
for such receivers with a noise temperature of 85° K. has been obtained as 0-37 W., 12:3 kW., 3-7 kW and 453 kW. for the 


Moon, Mars, Venus and Jupiter respectively. 


I. INTRODUCTION 


THE first successful attempt in the direction of inter- 
planetary radio communication was reported by 
Mofenson'. He used a radar set of 3 kW. peak power, 
transmitting 4-sec. pulse at 111-5 Mc./sec. and received 
the echoes from the Moon. Similar results have also 
been reported by De Witt and Stodola®. Clarke* 
discussed the possibilities of Earth-to-spaceship, planet- 
to-planet, and spaceship-to-spaceship radio communica- 
tion. Kerrt made a theoretical speculation on the 
possibility of getting echoes from the Sun and other 
planets at 30 Mc./sec. He has indicated the possible 
use of higher frequencies in the case of planets. The 
work reported here has been done with a view to finding 
suitable frequencies for such communication between 


different planets. Our main interest lies in using planets | 


as one-way radio communication circuits or as relaying 
reflectors rather than getting echoes. 


Il. GENERAL CONSIDERATIONS 


A signal can only be detected when the power reaching 
the receiver is such that its ratio to the Johnson noise is 
at least equal to the receiver noise figure. The received 
signal power at a particular point will depend upon its 
distance from the transmitter, the transmitting and re- 
ceiving antenna gains, and the molecular absorption of 
electromagnetic energy by gases in the path. 

Total absorption of radio waves (microwaves) in 
various planetary atmospheres has been reported 
elsewhere®. Power reaching a particular point can be 
calculated in terms of receiving and transmitting antenna 
gains in the following manner. 

The density of the power at a distance of r metres due 
to a transmitter which delivers a peak power of P, watts 
to an isotropic radiator is given by the total power 
divided by the surface area of the sphere of radius r. 

_However, if the transmitter is fitted with a directive 
antenna with maximum power gain G,, then the power 
density P, in the direction of maximum gain is: 





* Manuscript received 15 December, 1958, and in revised form 
12 November, 1959. 


JOURNAL OF THE BRITISH INTERPLANETARY SOCIETY 


G;. P, 
4nr* 


This power density is gathered by the effective absorp- 
tion cross-section of the receiving antenna, A m.*, and 
generates at the receiver input a power of magnitude 


P, = A.P,W. Xe eA (2) 


The absorption cross-section is related to the maximum 
receiver antenna gain G, by the following: 


G,.% 
2 
4g ™ ‘i me (3) 
where A is the wavelength of the wave radiated in metres. 


From (1), (2) and (3) we get the power at the receiver 
input as 


P,= W./m.? . * (1) 





— G,.Gy.P i 
' (4nrp 
The power of the transmitter which produces a power 
P., at the receiver is given by the expression: 
re P,.(4ary 
~ G,.G,% 
Now if the molecular absorption in the path is « db., 
then the power of the transmitter P’, required to produce 
the same power P, at the receiver is — 
P’, = P,.Antilog («/10) 


P,.A(4ar : 
= ae Antilog («/10) W. (6) 
By substituting the suitable values of G,, G, and 
minimum detectable power P,, we can have an idea of 
the transmitter power required for communication 
between Earth and various heavenly bodies. 
A different method of calculating the power of the 
transmitter is from the formula for the communication 
capacity of a channel*® 


S, 
C = W.log, (1 + ny) a (7) 


(4) 








P, . See 
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where S, = average received signals, 
W = channel bandwidth, 
Ng = noise power density, 
C = maximum or ideal rate of transmission in 
bits/sec. (if W is in c./sec.). 
Another useful relation which can be derived from (7) 
is: 
1-4 S, 
No 
Swerling*® has used Equation (8) to estimate the power 
requirements for space communication purposes. 


C< 





(8) 


Ill. CHOICE OF P, 
The minimum signal power which can be detected 
by a receiver can readily be calculated from the formula 


for noise figure : 
S/kTB 





NF = 9 
SulNo i" 

where NF = noise figure of the receiver, 
S = available signal power from signal source, 


S, = available signal power from receiver, 
9 = available noise power from receiver, 

Absolute temperature, 

Noise bandwidth of the receiver, 

= Boltzmann constant. 


ie 
| 


The minimum signal power which can be received 
corresponds to the minimum value of S,/No, which is 
unity. Substituting S,/N, = 1 in (9), the expression for 
NF reduces to: 


AY (P,)min- 
NF = ETB 


If the noise figure is given in decibels, 


(P, erin. = KTB. Antilog “*) (10) 


An average value of 15 db. has been taken for the 
receiver noise figure in calculations made here. 

Recently Townes ef al.’ have described the design and 
operating characteristics of a maser radiometer for use 
in radio astronomy at X band frequencies. The opera- 
ting system described has a bandwidth of 5-5 Mc./sec. 
and input noise temperature, including background 
radiation into the antenna, of about 85°K. The 
limitation on noise in such a system lies in the thermal 
radiation from components preceding the amplifier, in 
particular from the ferrite switch, the circulator and the 
hot ground seen by the antenna. A careful design may 
largely eliminate these noises, so that a noise temperature 
near 10° K. should be obtained. A noise temperature 
of 85° K. has been taken in our calculations. The noise 
figure is related with the input noise temperature by the 
following relation 


Input noise temperature = (NF—2) x 290° K. (11) 


In a practical case it is advisable to choose the value of 
P, = 10 X (P,)min. This is done to take into account 
the following factors : 


(1) Temperature used for the calculation of the 
Johnson noise is the average temperature (7,,,./+/2) 
of the planet and not the maximum temperature re- 
corded. Hence the instantaneous noise will be greater 
than the theoretical average value. 


(2) Noise level will increase when the planets come in 
line with a radiating body, such as the Sun. 


(3) Fluctuations and disturbances in the atmosphere 
due to sudden galactic bursts will cause reduction in 
signal intensity. 


It may be pointed out here that such a signal strength 
will be well above the total noise level obtained as a com- 
bined effect of Johnson, galactic and solar noises. 


IV. CHOICE OF ANTENNA GAIN 


Since the main interest is interplanetary communica- 
tion, it is important that the antenna structure be of such 
a type that long-period contacts are established. The 
situation allows only two types of structures, viz. (1) 
isotropic radiators, and (2) directive radiators which 
can be easily rotated to maintain the contact once 
established. The antenna gain should be as high as 
possible consistent with the limitation imposed by the 
second requirement. 

At 100, 500, and 1000 Mc./sec., high-gain arrays will 
be too big in structure and at a particular setting the 
transmitter will be in contact with the planet for a 
limited interval of time. The structure will have to be 
rotated, in order to maintain the contact with the planet. 
This is very inconvenient when the structures are very 
big. At these frequencies a simple isotropic transmitter 
antenna has been chosen (antenna gain 1). This has 
the disadvantage that a large proportion of the power 
will be wasted. The receiving antenna gains have been 
taken as 250, 500 and 1000 at 100, 500 and 1000 Mc./sec. 
respectively, since these can be realized in practice. 

In the case of 3000, 10,000 and 30,000 Mc./sec., the 
antenna structures will be more compact and rotation for 
continuous contact is no longer a problem. At these 
frequencies the antenna gain can reasonably be taken as 
2000. Owing to their high directivity these antennae 
avoid stray disturbances caused by delayed signals 
arriving as echoes from other planets. 


Vv. CHOICE OF BANDWIDTH 


Any value of the bandwidth can be taken depending 
upon the purpose for which the communication is made. 
The narrower the bandwidth, the smaller will be the 
power requirements for the transmitter (Equation (7)). 
In the present calculations a tentative value of 10 kc./sec. 
has been taken; this is a good approximation for speech- 
modulated transmission. 
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The bandwidth of 10 kc./sec. is required for speech 
transmission. It is most unlikely that this type of 
communication will be of a personal nature, making 
voice-to-voice contact a necessity. The teletype system 
of transmission is quite suited for the transference of 
technical data*. Above all, common experience of radar 
transmission shows that the signalling factor with pulsed 
power depends upon the peak power delivered and not 
the average power. This simplifies the problem of 
getting high average powers. 

In a very narrow bandwidth receiver, the setting of 
the receiver should take into account the Doppler shift. 
The Doppler frequency change is approximately given by 


Af=fV/V 


where V, is the relative velocity of the target and V is the 
velocity of the electromagnetic wave. 

Assuming the relative velocity of the Moon with respect 
to the Earth, due to the rotation of the Earth and the 
motion of the Moon in its orbit, to be 1000 miles/hr. an 
idea of the Doppler shift can be had from Table I. 


TaBLeE I.—Doppler Frequency Shift for Different Fre- 
quencies (Earth-Moon Communication) 








Change in frequency, 
Frequency, f, Mc./sec. A f, in kc./sec. 
100 0-15 
500 0-75 
1000 1-50 
3000 4-50 
10,000 15-0 
30,000 45-0 








VI. ABSORPTION IN THE PLANETARY 
ATMOSPHERES 
The data for the molecular absorption have been taken 
from an earlier paper by Chandra and Srivastava’ on this 


subject. Table II shows the absorption in the planetary 
atmospheres of different planets at various frequencies. 


TaBLE II.—Total Absorption of Radio Waves in 
Planetary Atmospheres 
































Absorption, db. 

Frequency, 

Mc./sec. Earth Moon Mars Venus | Jupiter 
400 0-26 ws me ah a 
500 0-01 ma sf Tg er 

1000 0-003 | — a ay pe 
~~ 3000~=—«| «0-092 | — | 0-002 | 0-003 | 0-505 
10,000 0-099 | — | 0019 | 0037 | 6374 
30,000 «=| 0160 | — | 0-087 | 0-170 | 29-498 

















Vil. RESULTS 


The data presented above have been utilised for 
computing the minimum power of a transmitter on Earth 
whose signals can be detected on various heavenly bodies. 
The results for an ordinary receiver are given in Table III. 
The results for a radiometer using masers are given in 
Table IV. 


TABLE III.—Power Required for One-Way Interplane- 
tary Communication Using Conventional Receivers 





Minimum power for the Earth’s transmitter 
(kW.) for one-way communication to 








Frequency, 

Mc./sec. Moon | Mars Venus Jupiter 
100 1:3 x 107% 43 x 10° | 1:25 x 10°| 1-35 x 16 
500 1-53 5-0 x 10° | 1-45 x 10° | 1-6 x 108 
1000 | 3-04 | 10 x 10 | 2-95 x 10° 3-5 x 108 
3000 70 x 10°; 23 x 10 7:0 x 10 | 85 x 10° 

10,000 8-0 x 10-77 26x 10° 7-5 x 10°) 36 x 108 

30,000 72 x 107?| 25 x 10° | 7-5 x 10°| 70 x 108 





TABLE IV.—Power Required for One-Way Radio Com- 
munication Using Receiver with Maser Preamplifier 





Minimum power required for the Earth’s trans- 
mitter (kW.) for one-way radio communication to 





Frequency, 


Mc./sec. Jupiter 


Moon Mars Venus 





x x | 
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10? 
10 
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S2e2 





* No maser operating at these frequencies has been announced. 


VIII. CONCLUSIONS 


The results presented in Table III indicate that radio 
communication with Jupiter under the conditions 
assumed in this paper requires a very large power which 
seems practically impossible. The higher microwave 
frequencies are attenuated owing to the presence of 
ammonia in its atmosphere. The power requirement may 
come within practical approach at a frequency of 3000 
Mc./sec. if the bandwidth of the receiver is further 
reduced. The power requirement for 100 Mc./sec. may 
also be brought within practical approach if a directive 
transmitting antenna similar to that at the receiver is 
used instead of the isotropic type, and if the bandwidth 
is further reduced. This will involve a very heavy 
structure which cannot be conveniently rotated for 
continuous contact. Thus communication with Jupiter 
at 100 Mc./sec., if made possible, will be only for short 
intervals. 

The power requirement of a transmitter for com- 
munication with Mars and Venus is least at 3000 Mc./sec., 
as compared to the power at other frequencies. Hence 
it is concluded that a frequency near 3000 Mc./sec., will 
be the most suitable for radio communication with these 
two planets. However, it may be mentioned that the 
power requirement at 100 Mc./sec. can be less than the 
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value at 3000 Mc./sec. provided an antenna with an 
approximate gain of twenty or more is used. Such an 
antenna will also have a fairly big structure and cannot 
be so conveniently rotated as a paraboloid at 3000 
Mc./sec. Further, the beam width of a 100-Mc./sec. 
antenna will be much larger than the angular diameters 
of the planets and thus some of the drawbacks of the 
isotropic antenna (disturbances due to delayed signals, 
etc.) will arise. 

For the Moon the power requirements are not very 
heavy at any of the frequencies and therefore any 
frequency may be used for communication. However, 
the use of 3000 Mc./sec. is recommended for the reasons 
already discussed above. 

Using a maser preamplifier, the power requirements 
will be much less because it can detect weak signals in 
noise. A comparison can be had from Tables III and 
IV. The calculations for 100 and 30,000 Mc./sec. have 
not been done because no masers operating at these 
frequencies have been announced. 

In all these discussions frequencies below 100 Mc./sec. 
have not been considered because at these frequencies 
the reflection and refraction effects in the ionosphere 
cannot be overlooked and it is probable that the signals 
may not reach other planets, as pointed out by De Witt 
and Stodola’. 


NOTE ADDED IN PROOF 

Since the acceptance of the above paper a large volume 
of information has been published in the April, 1960, 
issue of Proceedings of the Institute of Radio Engineers 
on Space Electronics. It has, therefore, become 
necessary to give a brief comment on the validity of my 
conclusions in the light of the information contained 
there. It may be stated that all the papers on this 
subject in this issue indicate that the best frequency 
region for space communication lies between 1000 and 
10,000 Mc./sec. and the frequency suggested in the 
present paper is 3000 Mc./sec. The following summary 
of the I.R.E. papers will be of interest in connection with 
the present conclusions : 


Mueller,® on the basis of an empirical relationship connecting 
mass, power, gain, and frequency of the parabolic antenna and 
transmitting system for space-vehicle-to-Earth communication, 
has concluded that the higher the frequency the greater will be 
the effective transmitted power. He takes the lowest possible 
frequency as 400 Mc./sec. and sets up an upper limit at 10,000 
Mc./sec., owing to increased problems associated with attitude 
control as the antenna bandwidth narrows and the increased 
mechanical tolerance at high frequencies. 

The spectral characteristics of various galactic and planetary 
radio noise sources are discussed by Smith,’° and it is suggested 
that the future of long-range communication lies in the micro- 
wave region. 

Ryerson" favours the use of passive satellite communication in 
comparison with the active system. The operating spectrum is 
assumed to lie between 2000 and 10,000 Mc./sec. The lower 
limit is chosen on the basis of the antenna gain and size as well 
as the complication arising from. possible diffraction from the 
passive system and the ionospheric effects. The upper limit is 
set by the accuracy of fabrication necessary to avoid phase 
cancellation and multipath effects as well as atmospheric attenua- 
tion considerations. 

The possibility of world-wide use of the “distributed” passive 


relay system is discussed by Bond ef ai.” in their paper “‘Inter- 
ference and Channel Allocation Problems Associated with Orbit- 
ing Satellite Relays.” The feasible operating spectrum is 
assumed to be 2000-10,000 Mc./sec. 

Dimond™ takes 1000-10,000 Mc./sec. as the available trans- 
mission “‘window,” on noise and atmospheric absorption con- 
siderations. The state of the art in present receiver and trans- 
mitter design has led to S band telemetry systems being favoured. 
A 2250-Mc./sec. telemetry system which is being attempted is 
considered. 

Brockman et al.* describe the TRAC(E) system for Pioneer III 
and JV, operating at 960-5 Mc./sec., chosen on the basis of solar 
and galactic noise, atmospheric absorption and the present state 
of the art. Calculation of the absorption of microwaves due to 
water vapour in the Earth’s atmosphere, showed that it was 
negligible at the recommended frequency of 3000 Mc./sec. Its 
contribution was approximately 5% of the total at 10,000 Mc./sec. 
and 30% at 30,000 Mc./sec. This is a further indication of the 
advantage of using a frequency of about 3000 Mc./sec. 
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LIST OF SYMBOLS 


A ee absorption cross-section of receiving antenna, 
B oo "bandwidth of receiver, Mc./sec. 

Cc maximum rate of transmission, bits/sec. 

G, maximum receiver antenna gain. 

G, maximum transmitting antenna gain. 

No available noise power from receiver, W. 

N,’ noise power density. 

NF noise figure, db. 

P, power density, W./m.? 


| receiver input power, W. 

(P,)min. Minimum detectable input power, W. 

P,, PP, transmitter power, W. 

Ss available signal power from signal source, W. 


So available signal power from receiver, W 
S, average received signals. 

T absolute temperature, ° K. 

Tmax. Maximum temperature recorded, ° K. 

V electromagnetic wave velocity, m./sec. 
Ve relative velocity of target planet, m./sec. 
WwW channel bandwidth, c./sec. 

f frequency, Mc./sec. 


Af Doppler frequency change, kc./sec. 
k Boltzmann constant. 

r distance, m. 

o molecular absorption, db. 

A wavelength, m. 


REFERENCES 

Mofenson, Electronics, 1946, 19, 92. 

H. De Witt and E. K. Stodola, Proc. Inst. Radio Engrs, 
N.Y., 1949, 37, 229. 

C. Clarke, J.B.J.S., 1948, 7, 49. 

J. Kerr, Proc. Inst. Radio Engrs, N.Y., 1952, 40, 660. 
Chandra and B. P. Srivastava, Z. Astrophys., 1959, 47, 
Swerling, Inst. Radio Engrs, N.Y., Trans. Mil. Electronics, 
1958, 2, 20. 

_ A. Giordmaine, L E. Alsop, C. H. Mayer and C. H. 
Townes, Proc. Inst. Radio Engrs, N.Y., 1959, 47, 1062. 
. O. Smith, J.B.1.S., 1953, 12, 13. 

.E. Mueller, Proc. Inst. Radio Engrs, N.Y., 1960, 48, 557. 
. G. Smith, ibid., 593. 
bs Ryerson, ibid., 613. 

. E. Bond, C. R. Cahn and H. F. Meyer, ibid., 608. 

.H. Dimond, ibid., 679. 

. H. Brockman, H. R. Buchanan, R. L. Choate and R. L. 
Malling, ibid., 643. 


10. 


pital a oP _ Yes th 


© The British Interplanetary Society. 1960. 








RELATIVE EFFECTS OF PROPELLENT DENSITY AND OF SPECIFIC 
IMPULSE ON ROCKET PERFORMANCE* 


By IRVIN G. HENRY,*+ M.S., Fellow 
(Communication from Aerojet-General Corporation) 


ABSTRACT 
The relative importance of propellent density and of propellent specific impulse is analysed for three different design 


situations. 


In each case a simple relationship for best performance is given in terms of the rocket component masses. The 


best mixture ratio of a liquid propellent (or formulation of a solid propellent) is determined. 


I. INTRODUCTION 


THE most important performance parameters of a 
propellent are its specific impulse /,, and its density 
p. If one propellent is superior to another in both of 
these respects, it will certainly yield a better missile 
performance. However, if the density is greater and 
the specific impulse less, or vice versa, the effect on the 
performance is more problematical. It has become 
customary to select propellents for maximum p"/,,, 
where n may vary from 0 to about 1-5. The following 
analysis will show how the correct value of m is to be 
chosen. 

Three design situations will be considered; in all of 
them the rocket will accelerate rapidly to a high burnout 
speed and will then coast ballistically or aerodynamically ; 
cruise-type rockets are therefore excluded. In the first 
situation, part of the mass of the empty rocket is fixed 
and the propellent volume is fixed. In the second 
situation part of the mass of the empty rocket is fixed 
and the full mass is fixed. The third situation allows a 
varying diameter, but fixes the full mass, the length and 
part of the empty mass. 

In the first two cases the intention is to maximize the 
burnout speed v of the rocket, which is given in terms 
of the propellent jet speed (j) and the masses of the rocket 
when full and when empty (m, and m,, respectively) to a 
first approximation by: 


v = j In (m,/m,) Ga . (1) 


The application in these two cases can be to a ballistic 
missile, or to an air-to-air boost-glide missile with fixed 
diameter ; in both cases the range is an increasing function 
of boost speed, so it is sufficient to maximize the latter. 
The third case involves a boost-glide missile of variable 
diameter. Here also the intention is to maximize the 
glide range, which is roughly proportional to the product 
of the empty mass, some power of the boost speed, and 
the inverse square of the missile diameter. 


Il. PROPELLENT VOLUME AND PART 
OF EMPTY MASS FIXED 


The first case we shall consider is a rocket with an 
empty mass composed of two parts. The first part, 
Mx, is fixed ; this might correspond to the payload and a 
portion of the propulsion system. The second part, 
corresponding to the propellent chamber, has a mass 
M,. proportional to the propellent volume V. The 
volume V will be taken fixed, so that m,, is fixed, and 
therefore also m,. The propellent has density p, and 
mass m, = pV. Therefore 

3 My om, + pV 
v=Jln a jin ~% 


(2) 


If the propellent undergoes small changes dp in density 
and dj in jet speed, the corresponding change in burnout 
speed dv is given by the relation 
m+pVi... mM, V 
dy = | In ——— | dj + j ————— ' — 
v [ “ orig me (3) 
If the effects of dp and dj cancel, so that there is no 
change in v, we have 
m, Mm, 
dj m, m, 4d 
1. mm 5 . 
ji In(m,/m,)_ p 
We must now find the value of n in 
p"!,, = constant ss ihe (5) 
which gives the same result. We rewrite (5) 
p"j = constant 
and differentiate the logarithm of this expression, 
obtaining 


dj d 
Fi nta oe noel 
J p 
Comparing (4) and (6), it is clear that 
oes. * een me 
" In(eny]m,) a 
In{ — 
1—x 
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where x is m,/m;. It is remarkable that the value of n 
should thus depend on the vehicle propellent fraction 
alone, and not directly on the density of the propellent. 
Selected values of n(x) are given in Table I. 


TABLE I 





x vee ae AS 
4 In( i ) 
1-—x 
1-00 
0-898 
0-783 
0-721 
0-655 
0-497 
0-390 
0-317 
0-0 
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A hypothetical example of a rocket designed under the 
constraints given might be an air-to-surface missile of 
fixed payload mass and payload diameter, weighing 
1000 Ib., of which 200 Ib. is propellent. It is desired to 
maximize range, but instead speed is maximized, since 
range is an increasing function of speed. Now x = 0-2, 
therefore n=0-9. Accordingly it is mecessary to 
maximize the quantity p®-*/sp. Consider next a single- 
stage ballistic missile with a specified payload designed 
to fit within a predetermined volume, such as a tube in 
the ground or in a ship. A fixed volume remains for 
propellent, so that it corresponds to the constraints given 
in this section. Such missiles normally have larger 
propellent fractions and smaller values of nm. For 
example, if the missile weighed 10,000 Ib., of which 
7000 Ib. was propellent, it would be necessary to 
maximize p®-**/,,.. 


Ill. FULL MASS AND PART OF EMPTY 
MASS FIXED 


The second case to be analysed is a rocket in which, 
as before, the empty mass comprises a fixed part m,, 
and a part m,. proportional to the propellent volume V. 
However, the full mass m, will be fixed, while the 
propellent volume will not. Thus 


Mee = (k/p) mM, = (k/p)(my a Mea) ‘* (8) 
or 
_ (k/pmy, — Ma) 
Meo = “th ip) ‘ - (9) 


In modern design the ratio of chamber mass to propellent 
mass, k/p, is rarely greater than 0-15. Accordingly we 
may neglect k/p in the denominator of Equation (9). 
Therefore, differentiating (9), one finds 


dm,s =(- k/p*)(m, "1 Ma) dp 
=—Madpip .. 3 ve - 


but 
v = jin (m,/m,) 


and since dm; = 0, we have 





a . m, _. dm, 
d= din + (11) 
Because dv = 0, we have also 
di 
J on Cis .. aD 


ij a m, In (m;/m,) zs 


But m, = my + m,, and dm, = dm,,. Substitution of 
(10) into (12) therefore gives: 
dj dp * Mes 


> A s.. Seam oe 
j i p m, In (m;/m,) (15) 


Therefore, in this case 


Me2/M, 
a= ingypin ae .. (14) 
A typical single-stage ballistic missile might have m,,./m, 
=0-5 and m,/m,=5. This would yield nm = 0-31. 
Thus, as is to be expected, when the missile is not 
volume-limited, the value of n is smaller; i.e., the burnout 
speed depends less on the density p. 


IV. FULL MASS AND PART OF EMPTY 
MASS FIXED, VARIABLE DIAMETER 


In the third situation to be analysed, we wish to 
maximize the glide range of a missile. As before, part 
of the empty mass is assumed fixed; this might corres- 
pond to the payload. The full mass is also fixed. 
However, the third problem will differ from the second 
one considered in that the missile diameter is variable 
and the length fixed. 

The glide range s is assumed proportional to the empty 
mass, a power of the burnout speed and inversely propor- 
tional to the frontal area; that is 


y’ 
s| og im jaey{ » — wy} «vy fe 


v 


For pointed missiles it is found that p is about 0-6 
since drag varies as y®~°-* — y'-*; for blunt-nosed missiles 
p is about 0-1 or less. From missile geometry we have: 


d? x (my/plmKImily) -.  -» (16) 


where /,, is the missile length and /, is the length of the 
propellent chamber. It was specified that /,, was 
constant. The fundamental rocket equation gives 


v = jin (m,/m,) ob in (2) 


We now make the approximations that /,,//, is constant 
and that (v’)” is much less than y”; the second approxi- 
mation can be roughly satisfied by adjusting the value 
of p. These approximations, together with equations 
(2), (15) and (16), yield 


? 
s & (m,/m,) pj? [in (mim, | -» (17) 
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We write 
mM, = My — mM, = My — Ma — My -. (18) 
and 
dm,=dm.. .. es .. (19) 
Thus 


a(™) = — = a cb Ce 


mM, m, m, m, m, 
or, on using Equation (10) 
a(~) aS ae 
m, m me p 


Rather than maximize s, we prefer to maximize s''”. 
Doing this, and substituting Equations (18) and (21), 
we find 


i My Meg Mes 
din si /P) = | 1 bs, m, Mm, m, dp ) dj 
3 af ne |r J 
a me m,.— 
=0.. a ae = “43 "aa 


Thus finally 
— l {1 a. cima) | i M,2/M, (23) 


Pp (m,/m,) — 1 In (m;/m,) 


In a typical case one might have m,/m, = 2 (i.e., pro- 
pellent mass is half the missile mass), m,. = 0°25 (i.e., 
propellent chamber is one-quarter of the empty missile 
mass), and p = 0-6 (missile with pointed nose). We find 
n= 1-19. This is a large value for n, reflecting the 
importance of reducing the frontal area. 





Vv. SUMMARY AND CONCLUSIONS 


We have presented a method for selecting the best 
propellent or the best mixture ratio of a propellent 
combination for each of three design situations. It was 
assumed necessary to select the propellent having the 
maximum value of p"/,,, where m is appropriately 
calculated. The expressions for m in the three cases 
considered are given in Equations (7), (14) and (23). 

When nv has been computed, it is convenient to use it 
in conjunction with a plot of the type shown in Fig. 1, 2 
and 3. Here the logarithm of /,,, is plotted against the 
logarithm of the density, so that contours of p"/,, 
appear as straight lines, with a parallel family for each 
value of n. Each bipropellent combination generates 
an arc, on which the best mixture ratio is that which 
gives the point furthest up and to the right on the grid 


having the applicable value of m. Monopropellents 
and solid propellents generate points rather than arcs, 
but the choice is made on the same basis. 

The use of the Figures is easily illustrated by example. 
Consider a single-stage ballistic missile meeting the 
conditions of Section III with m,/m, = 10 and m,./m, = 
0-5. This gives nm = 0-22. Consulting Fig. 1, we see 
that, with this slope and using oxygen/hydrogen, the 
optimum mixture ratio (oxidizer/fuel) is 5-7. It is also 
apparent that in this case fluorine/ammonia gives con- 
siderably better performance than oxygen/hydrogen, 
because of the higher density of the fluorine/ammonia. 
Since these numbers are typical of modern rockets, one 
must conclude that fluorine/ammonia has been excluded 
on grounds of logistics rather than performance. 

Some practical aspects of the foregoing derivations 
should be pointed out. First, the application to a multi- 
stage rocket requires an extension of the theory just 
given; for example, a change of propellent in the second 
stage might increase the speed of that stage alone but 
also it might reduce the speed of the first stage. Secondly, 
the best choice of propellent is, at times, quite different 
after the vehicle design is frozen from what it was before. 
This frequently occurs when a rocket is designed for 
fixed gross mass; then later the restriction on gross mass 
is lifted but the motor case dimensions are retained, thus 
imposing a restriction on volume. Thirdly, we should 
warn against any oversimplified comparison of solid and 
liquid propellents because (a) the effective density of 
solids must be reduced by 15 to 30% to allow for the 
presence of the port, (6) each type of propellent affects 
the empty mass in a different way ; a solid rocket requires 
a larger safety factor in chamber wall thickness and may 
leave more unburned propellent. Therefore solids and 
liquids should not be compared directly on an J/,,, versus 
p diagram. 


LIST OF SYMBOLS 


I,» specific impulse of propellent. 

j jet speed of propellent. 

k ratio of hardware mass to propellent volume, m,,/V. 

1, missile length. 

1, propellent chamber length. 

m, empty rocket mass. 

m. part of empty rocket mass proportional to propellent 
volume. 

Me fixed part of empty rocket mass. 

m, full mass of rocket. 

m, propellent mass. 

n exponent of propellent density. 

Pp speed exponent in equation for glide range. 

s glide range. 

v burnout speed of rocket. 

v’ speed of rocket at end of glide. 

V_ volume of propellent chamber. 

p __ propellent density. 
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AN ENVIRONMENTAL CONDITIONING SYSTEM FOR A MANNNED SATELLITE* 


By NORMAN J. BOWMAN,+ Ph.D., Fellow and EDWARD H. DINGMAN? 
(Communication from The Martin Company) 


ABSTRACT 


A specific design is presented for an environmental conditioning system for use in a manned satellite designed to stay 
aloft for 60 days. Water is recovered by evaporative distillation and refrigeration, oxygen regenerated by electrolysis of 
water, carbon dioxide removed by use of lithium oxide (non-regenerative), and capacity for radiating 150,000 B.T.U./day 
from the satellite as a whole provided. With automatic control this system is suitable for use in an animal-carrying satellite. 
The weight is favourable when compared with other systems that have been suggested, and well within the conditioning allow- 


ance of a 5000-lb. manned satellite. 


I. INTRODUCTION 


THE problem of controlling the environment in a manned 
satellite, which is to remain aloft for more than a few 
days, is a difficult but essential one. It is necessary to 
remove heat, control humidity, remove carbon dioxide 
and replenish oxygen. Many relatively simple ways of 
doing this are known, but they universally suffer from the 
defect of high weight and power requirements. 

A hypothetical case has been chosen to permit the 
making of calculations for a system of control which is 
applicable to a wide range of purposes. We have chosen 
the case of a satellite with an approximate 90-minute 
orbit carrying a man and designed to stay aloft for 60 
days before recovery. Such a satellite would not need to 
weigh over 5000 Ib. and is within the limits of accom- 
plishment in the near future by use of the Titan missile 
as the basic booster. Power will be supplied to the 
satellite by solar cells with an average daylight level of 
approximately 14 kW. A sun-seeker is to be used for 
solar cell orientation and part of their output power is 
stored in batteries for use while the satellite is in shadow. 
The power from the solar cells must satisfy requirements 
for all purposes including control, telemetering and 
environmental conditioning. Peak loads, as for example 
communications and command data telemetering, are 
satisfied with the assistance of stored battery power. To 
meet requirements, an estimated 120 ft.? of cell surface 
area will be needed. 

Heat removal is one of the more serious problems in a 
satellite, as not only the cabin but the electronics com- 
partment must be cooled and the heat radiated to space. 
The design presented here is capable of handling 150,000 
B.T.U./24 hr. in the life-ambient range of temperature. 

The fact that the satellite is manned simplifies the 
control system and increases general reliability. The 
same system could be used for an animal-carrying bio- 
logical satellite with the addition of a number of 
automatic control devices. 

Details on the major parts of the system are given 
below. 


II. ATMOSPHERE CONTROL SYSTEM 
1. DISCUSSION 


The atmosphere used in this system will consist of 
144 mm. of oxygen pressure, 130 mm. of helium pressure, 
and not over 5 mm. of carbon dioxide pressure. The 
water vapour pressure is variable. Table I gives the 
quantity of various products a man will produce or 
need'~* under light activity in 24 hours. 








TABLE I 

Product or | Rate of Production | 

Requirement | or Use Remarks 
Water 6°6 Ib./day | Produced (less required) 
Oxygen .. oad 2-2 Ib./day Required 
Carbon Dioxide | 2-2 Ib./day | Produced 
Methane oe 0-3 1./day | Produced 
Heat (biological) | 12,000 B.T.U./day | Produced 
Heat... ..| 150,000 B.T.U./day | Total to be radiated 


from satellite including 
all sources. 





The problem of the oxygen supply was examined 
quantitatively. The use of gaseous oxygen leads to a 
very high mass requirement (about 300 Ib. for the 
oxygen and tanks alone) for a 60-day satellite, largely due 
to the mass of the tanks. Utilization of liquid oxygen is 
quite feasible using lightweight vacuum-insulated con- 
tainers at about 2 atm. pressure. The evaporation rate 
under such conditions should not exceed the rate at which 
oxygen is needed initially, and with the best insulation 
external heating would be needed to maintain the neces- 
sary evaporation rate. A mass of 150 Ib. total is 
estimated for the oxygen supply if this system is used. 

It is apparent that a man will produce, because of 
combustion of foods, more water than he consumes. It 
is desirable to reclaim essentially all water available. 
This is discussed in a later section. Once recovered, 
this water over and above that needed for drinking 
purposes is surplus and available for decomposition to 
recover the oxygen. This would result in a substantial 
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weight saving in the atmosphere system. The only 
oxygen that need be lost to the system is that going into 
carbon dioxide. This amount minus the amount pro- 
duced by food combustion is about 60 lb. We need, 
therefore, to carry about 75 Ib. of water at takeoff in 
order to have a balanced system for 60 days and an 
adequate safety factor. If waste water is not reclaimed, 
300 Ib. of water in addition to 150 lb. of oxygen plus 
container would have to be carried. 

A system has been devised whereby the oxygen needed 
is produced by electrolysis of water. It utilizes a low- 
speed rotating cell where the “centrifugal force’ keeps 
the gas and water separate. 

In order to conserve battery power, the electrolysis cell 
is run only during the time the satellite is in sunlight and 
direct power from the solar cells is available. Enough 
oxygen is produced and stored during this time to last 
while the satellite is in the Earth’s shadow. The hydro- 
gen formed in the electrolysis can be vented to space or 
used to replenish the cooling system gas if it is used as 
heat transfer medium. 

Another problem is that of keeping the carbon dioxide 
content of the atmosphere below 5 mm. partial pressure. 
For a trip of this duration and with only solar power 
available, the use of a regenerative system, such as that 
based on calcium oxide, is not economical. Lithium 
oxide is the best non-regenerative agent available from 
the standpoint of mass. To remove the carbon dioxide 
produced, a total of 80 lb. of lithium oxide will be 
required. A 25% safety factor has been added, making 
a total of 100 Ib. 

Heat removal from the cabin or animal cells proves to 
be the controlling factor from the standpoint of air 


re-circulation. In order to remove the 12,000 B.T.U./24 
hr. produced by a man or animals, an air circulation 
rate of 80 ft.*/min. is required if air is to enter at 58° F. 
(14-4°C.) and leave at 80°F. (26-7°C.). This would be 
about 80% of available air exchanged each minute if a 
compartment of 100 ft.* is used. Since the removal of 
carbon dioxide and water does not need to proceed at 
such a high rate (with respect to air flow), an economy 
in power and weight is effected through by-passing 
70% of the air flow around the carbon dioxide and water 
removal system. 

The liberation of body gases presents a problem. 
The major constituent is methane and in 30-40 days a 
concentration would build up that could reach the lower 
explosive concentration of 5% by volume of the cabin 
atmosphere. It is therefore necessary to remove the 
methane. Unfortunately this is not simple as it is an 
extremely stable molecule. One possible way is to 
rapidly and completely exchange the entire atmosphere 
at the end of 30 days. However, by use of a small 
platinum coil heated electrically to 500°C. it is possible 
to burn the methane to carbon dioxide and water. This 
is carried out just before the gas enters the blower and 
after the oxygen content has been raised to 144 mm. 
partial pressure. Any hydrogen present would also be 
oxidized. As methane generation is slow and a con- 
centration of several percent is tolerable, it is not neces- 
sary to cycle the entire atmosphere stream through this 
step. The rate of air flow can be maintained through 
this cell at 2 ft.3/min. and the remainder of the atmosphere 
97-5%) by-passed. The power requirement and heat 
produced with this low flow rate in a cell using a heat 
exchanger is quite small, approximately 10-15 W. 


Temperature and humidity 
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Fic. 1. Atmosphere system. 


Pressure: 5 Ib./in.? Circulation : 24 ft.*/min. in Condenser, 80 ft.*/min. total. 
Humidity : 40-50%. Composition: O, 144 mm., He 130 mm., CO, 7 mm. (max.), H,O variable. 
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Other body gases (mainly sulphur compounds) are 
removed by a small activated charcoal bed. 

Oxygen and nitrogen free-radicals formed by space 
radiation and particle collisions with the cabin gases 
could prove a health problem if allowed to concentrate. 
However, products, mainly ozone, are destroyed both 
by the activated charcoal-silver bed and by the catalytic 
burner. It is likely that some of the ozone will react 
with the methane while still in the cabin volume. If 
helium is used instead of nitrogen, nitrogen free-radicals 
will not be of concern. This might simplify environment 
control. 


2. DESCRIPTION OF THE SYSTEM AS SHOWN IN Fic. I. 
A flow diagram of the specific system is shown in 
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mixes with the gas in the chamber and leaves by screened 
vents. The gas temperature will rise from 58°F. 
(14-4°C.) to 80°F. (26-7°C.) while passing through the 
compartment. The temperature, humidity, and oxygen 
content of the incoming air is automatically monitored 
and is subject to automatic or manual control at the 
operator's discretion. The 80°F. air forced out of the 
cabin passes through a heat exchanger (Fig. 2) which 
cools it to about 56°F. (13-3°C.). At this point, approx- 
imately 70% is by-passed and returned to the cabin 
without moisture or carbon dioxide removal. 

The remaining 30% of the gas at 56°F. is passed 
through another heat exchanger where it is cooled to 
36°F. (2:2°C.). This results in a mixture of saturated 
air and water droplets which is passed into a centrifugal 
separator (Fig. 3) where water is removed and sent to 



























Fig. 1. Gas of proper composition and temperature 
enters the compartment under flow rate control; It storage. A vortex separator using “moisture lock” 
Outlet manifold BAN tube 
Thin metal sheet with \_H, outlet fittings 4; inlet 
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Hydrogen flow am. ft. 3/min. ae 60 60 
Air flow rate, ft.*/min. .. 80 22 
Air velocity, ft./min. - wt 4 2 
Hydrogen inlet temperature, © F. 34 23 
Hydrogen outlet temperature, F. 52 34 
SECTION B-B Air inlet temperature,° F.  .. 80 56 
Air outlet temperature, ° F. 56 36 
Required heat removal by hydrogen, B.T.U. ‘/min. 9-4 5-0 
H,O gas from Diffuser plates,5 req'd 
heat exchanger ‘ A— 
AN type fittin ye deg 
Fic. 3. Water separation cell. Pt “a Slip joint 
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Air temperature, ° F. 36 36 
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filler may be used in place of the centrifuge. The 
saturated gas at 36°F. is then passed into the carbon 
dioxide removal chamber (Fig. 4). It flows first around 
the outside in a double wall to keep the absorption 
chamber cool and to pick up heat (to reduce the relative 
humidity). The amount of heat liberated from he 
overall reaction : 
Li,O + CO, = Li,CO, 

is sufficient to heat the normal flow of gas to about 
56°F. The lithium oxide is in granular form and is 
separated into 2-in. thick layers by perforated plastic 
disks to minimize caking and give relatively unimpeded 
gas flow at all stages of the tube’s life. The air will 
contain the moisture necessary for the above reaction to 
take place. The rate of the reaction: 

Li,O + H,O = 2LiOH 
is slow at medium temperature, and neither the carbonate 
or hydroxide is hydroscopic to any extent. This system 
will keep the carbon dioxide at a low level, certainly well 
below the 5 mm. maximum. 

The exit gas flows through a charcoal bed to remove 
other gases. From here it joins the by-pass line and the 
joint stream passes through the automatic oxygen 
analysis and make-up system. Most of the gas then 
passes through the blower, by-passing the catalyst 
chamber for methane removal. The 2-5% which does 
pass through this chamber will pass over a coil heated to 
500°C. and the methane and any hydrogen that has 
leaked into the system oxidized. Possibly the cold 
oxidation of methane and hydrogen can be brought 
about by a catalyst, such as divided platinum, combined 
with one or another radioactive isotope. An emer- 
gency oxygen and helium supply ties in just before the 
blower in the air cycling system. This supply is suffi- 
cient to allow an emergency descent of the vehicles to the 
surface. 

The main oxygen supply is furnished by the water 
electrolysis system and is shown in detail in Fig. 5. To 


supply the 2:2 lb./day oxygen which a man needs will 
require the electrolysis of about 2:5 Ib. water. In 
practice, about 2-86 lb. water per 24 hours will be 
electrolyzed. The electrolysis of water can be carried 
out at approximately 95° efficiency at moderate current 
densities according to date available‘. Approximately 
190 W. constant power is required. However, no solar 
power will be available when the satellite is in shadow; 
this is about one-third of the time for an equatorial orbit 
of near 90 minutes. To conserve battery capacity 
requirements, it is planned to stop electrolysis during 
dark periods; therefore about 285 W. for 15 hr. a day 
will be required. The actual cycle will consist of shut- 
down for approximately 30 min. and operation for 60 
min. out of every 90 min. The oxygen needed during the 
periods of shutdown will be stored at one atmosphere 
pressure in the low-mass shell surrounding the cell. 
Hydrogen liberated at one of the electrodes will either 
be vented to space or used in the cooling system. The 
electrolyte will probably be acid rather than the more 
commonly used alkali because of emulsion problems. 
This point, however,must be checked experimentally. 
The problem of rotary gas seals for this cell was solved 
several years back by refrigerant pump manufacturers. 
A loss of one part in 10* per year is considered excessive 
for such seals. However, it appears entirely feasible to 
do away with the rotating cell and use glass fibre for 
holding the electrolyte. Some research needs to be 
done on such cells relative to free-fall environment. 
One part of the cycle that will be operated manually 
and intermittently is the recovery of moisture from urine 
(1500 ml./day). Feces and urine will be placed in 
separate plastic containers and the urine returned into 
the system as shown in Fig. 1. The urine is introduced 
into the part of the system beyond the by-pass and before 
the exhaustive treatment given the residual 30% of the 
gas. The evaporator is fitted with a by-pass system 
which will be the lines used except during evaporation. 
Periodically the urine will be well filtered and then 
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B YS Inlet air space between 
two Plexiglass cylinders 


SECTION A-A 


| Fic. 4. Lithium oxide container. 
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Lithium oxide. . 100 Ib. (25% excess) 
Be Heat generation 0-8 kcal./min. 
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Fic. 5. Electrolysis cell. 


2-6 Ib. water/24 hr. 


190 W. constantly or 
285 W. for two-thirds 


Electrolysis rate 
Power requirements .. 


of the time 
System speed . ~10 rev./min. 
Sphere diameter 21 in. 


oe 





Programmer and 


sensing system 










CaO with platinum 
catalyst no scale 














| Fibreglass sphere 18 Hydrogen demand valve 
2 Electrolysis cell SECTION VIEW 19 Motor and drive gears 
3 Oxygen shaft 20 Bearing 

4 Water supply cube 1! Oxygen outlet tube 21 Rotating contacts 

S Electrolysis plates 12 Demand valve 22 Rotary seals 

6 Separator plates (bottom perforated) 13° Rotary seal 23 Hydrogen supply tube 
7 Separator plates (bottom perforated) 14 Bearing 24 Check valve 

8 Water and acid 1S Rotating contacts 25 Hydrogen shaft 

9 Stationary agitator _ 16 Pressure sensor 26 Oxygen demand vaive 
10 Water and acid level sensors 17 Pressure sensor 27 “Perforated support 


sprayed slowly into the evaporator. Urine bacteria 
can be controlled by a trace of silver salt. This evapora- 
tor consists of a series of perforated silver-plastic baffle- 
plates arranged so that considerable turbulance results. 
The droplets are broken up into fine spray which will 
evaporate rapidly leaving the filterable solid materials 
in the chamber. The silver-plastic baffle plates can be 
periodically replaced as solids build-up dictates. These 
plates will have a mass of some 4 oz. per set. As the 
feces normally contains little moisture (about 150 ml./24 
hr. is produced) it appears simpler to discard it than to 
remove the moisture. If necessary this moisture could 
be recovered but at some mechanical difficulty. The air 
in this part of the system will normally bypass the evapor- 
ator except for periods when urine is being processed. 
The moisture from the urine and/or feces will be a part 
of the 30% of the gas which is treated exhaustively to 
remove moisture, carbon dioxide, gaseous acid con- 
stituents, micro-organisms, and other impurities as 
indicated in the flow diagram. The result will be odour- 
less and sterile air and water of satisfactory purity. 
The water will be the equivalent of distilled water and 
suitable for drinking or electrolysis. 

A weight and power requirement breakdown is 
summarized in Table II. 


Ill. WATER CYCLE 


The overall cycle for water is shown in Fig. 6. The 
water is forced on demand by the gas pressurized dia- 
phragm into the cabin. The water will leave mainly as 
vapour or as urine. The recovery of urine and moisture 
in the air as liquid water has been discussed previously 


as it is integrally connected with the atmosphere cycle. 
After leaving the centrifugal separator the water passes 
through a silver-charcoal bed to remove gas and any 
materials carried over and to destroy bacteria. Next, 
the treated water goes to storage. The storage tank 
also holds replaceable activated charcoal inserts to fur- 
ther remove dissolved gas. The water produced is 
completely satisfactory for consumption although there 
is admittedly a psychological problem in so convincing 
a man. 





CABIN 















Gas and water vapour 
mixture (see Fig. 2) 


H,O storage 


Charcoal 


Gas pressurized 
and silver Separator 
bed 


diaphragm 


Heat exchanger 






To electrolysis cell 
Fic. 6. Water system. 


The passenger may drink as much or as little as he 
desires. In one way or the other the water is always 
recovered and the only oxygen lost to the system will be 
that tied up in carbon dioxide as previously mentioned. 
In 60 days this loss to carbon dioxide will equal approxi- 
mately 150 lb. water. About 75 Ib. water will be aboard 
at launch, the balance being supplied by the biological 
utilization of food. If the man’s food supply is not 
completely dehydrated, there is a further net gain or 
build up of water in the system which has not been con- 
sidered here, that weight of food which is trapped water 
can be deducted from the 75 Ib. of water carried. Water 
used by the occupant for washing purposes is similarly 
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TABLE I.—Conditioning System 

















Equipment Mass, Ib. Power, W. 
Atmosphere Control Equipment 

Urine storage plus ee filters 12 — 
Feces storage . 4 — 
Urine evaporator 5 _ 
Solid particle filter (including glass wool) . 4 _ 
Heat exchanger (air coo 12 (15 for He) — 
Heat exchanger (water condensation) 8 (12 for He) — 

Centrifugal separators - 5 Same shaft as 

electrolysis cell 
Battery to support conditions ey ++] 10 _— 
Carbon dioxide mete P ++ 110 — 
Charcoal bed 4 5 _ 

Electrolysis cell 8 *285 max. 

(zero min.) 
Initial charge of water in electrolysis cell 10 — 
Motor to run cell 1 25 
Combustion chamber (catalyzed) a 2 20 
Blower .. , ; *125 
Emergency Oxygen Tank plus ‘Oxygen gas —_ 
Emergency helium tank plus helium 10 (50 for He) — 
Miscellaneous lines, controls, regulators, ete. 20 

234 (281 for He) 475 max. 
380 average 


Water System 
Contact chamber, heat exchangers and centrifugal separators 


Charcoal-silver bed. bi 
Water storage chamber, etc. .. 
Water supply, initial .. 


Heat Radiator System 
Heat exchange gas blower 
Black body radiator 
Bladder valves w/control valves, lines, etc. 


10 Ib./in.? compressor and 1 ft.* tank eae 1 ft.2 of Bas per hour) s 


Fittings, lines, connectors, etc. 


Je Mass and power are entered under the 
atmosphere system. 


5 
10 | es 
.| 75 _ 
x | 
| | 
e 5 (8 for He) 265 (575 for He) 
eal 37 (50 for He) _- 
a 8 (20 for He) 1 
4 10 
50 poate 








276 (586 for He) 





Totals 


| 

b S 

| 104 (132 for He) | 

wo) 428 (503 for He) Max. 751 (1061 for He) 
| Average 656 (966 for He) 





* No power is required in shadow as cell is shut down to save battery weight. Time in shadow is assumed as being one-third of time. 


re-cycled and, within reason, he may use as much as he 
likes since the cooling system is more than adequate to 
recover it. The wash water would, of course, be 
re-cycled with the urine as a flush. As it is awkward 
under free fall conditions to handle fluids in the open, it 
is expected that this problem will be circumvented 
both by moving free fluids about in polyurethane 
sponges and plastic squeeze containers. Water droplets 
escaping will evaporate or be swept by the moving air 
into the exhaust outlets. 


IV. REFRIGERATION SYSTEM 


Assuming that over 1-5 kW. power is taken in from 
solar cells and that this plus body heat and solar leakage 
into the satellite must eventually be disposed of, it 
becomes a matter of necessity to consider the method of 
heat removal early in the system design. Only dumping 
of mass at high temperature or radiation cooling are 
available as methods for heat removal from the satellite. 


A quick calculation shows that mass-dumping is im- 
practical because of the quantity of mass (preferably 
water) which would have to be carried. 

A total of as much as 150,000 B.T.U./day (a consider- 
able safety factor is contained in this number) from all 
sources must therefore be removed by radiation. Calcu- 
lation shows that to radiate this amount of heat to space 
requires about 110 ft.2 of black body surface if the 
temperature differential from 5°K. in space to 353°K. 
(80°C.) in the electronic chamber could be utilized. A 
total of 142 ft.? of radiator is used to provide adequate 
safety margin. The design of the radiator itself is shown 
in Fig. 7 and 8 and the overall refrigeration system in 
Fig. 9. The Appendix gives the thermodynamic details.*.* 

The radiator consists of ten separate channels running 
in alternate spirals about part of the hull. These 
channels all separate where the gas leaves the manifold 
and passes through the bladder valves. This is done so 
that in case of meteor holing or leakage the channel 
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a 5°C. H, “seo ~ Fic. 7. Hydrogen cooling system radiator. 
f a joa NorTEs 
Radiator channels ~0°01 in. 
: Aluminium or stainless steel sheet. 


Channel width ~4 in. 
Channel depth ~1 in. 
Channel spacing ~0°5 in. 
Turns per channel ~1. 
Channel length ~30 ft. 
Channel volume ~1440 in.* 
Manifold rings ~96 in. 0.d., ~90 in. i.d., length 10 in. 
10. Bladder valve housing length ~3 in., diameter ~2 in. 
Out 11. Inlet diameter ~4 in., outlet diameter ~3 in. 
manifold 12. Gas flow per channel ~6 ft.*/min. 
13. Gas flow total ~60 ft.*/min. 
14. Inlet gas temperature ~80° C. 
15. Outlet gas temperature ~—5° C. 
16. Gas stay time ~8 sec. 
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Heat radiation to space 
H, in at 80°C. 2) H, out at -5°C, 


o7yra_Heat exchanger and 
C Blower ‘ ; rm H,O condenser 
H, supply from Pressure data ne ie 
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H,O electrolysis to telemetry r 


\cHeat exchanger 
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From hydrogen supply 


Circulation through for re-entry purge 


electronics compartment 
Purge znd pressure 
relief valve vented to space 


Fic. 9. Hydrogen system. 


Pressure $4 y 15 Ib./in.* 
Flow through radiator 60 ft.*/min. 
Radiator surface area 142 ft.? 


involved can be automatically shut off from the system 
before much gas has been lost. The system should be 
operable even if three of the ten channels are put out of 
commission. The radiator is made of two sheets of 
thin flexible metal brazed together to form gas channels 
which can be inflated by compressed gas. 

The electronics compartment is designed to run at 
80°C. exit gas temperature. Gas at ~ 15°C. is diffused 
in at one side and a blower forces out the heated gas at 
the other side. The entire system will operate at 760 mm. 
in contrast to the 275 mm. pressure in the cabin atmos- 
phere system. 

The gas at 80°C. will be cooled by radiation to —5°C., 
whereupon it passes through the two cabin heat ex- 
changers and then into the electronic and mechanical 
equipment compartment completing the cycle. As 
demand varies, depending on whether the satellite is in 
the sunlight or shadow, the coolant gas flow must be 
variable. 

The choice of a heat transfer gas was considered 
carefully and hydrogen chosen for several reasons. The 
first is its high heat capacity and low viscosity. The 
second is the fact that it can be replenished in case of 
leakage from the hydrogen produced in the electrolysis 
cell. A major disadvantage is the explosion hazard. 
The explosive limits are from about 5% to 94% by volume 
of hydrogen in oxygen. However, hydrogen has been 
used for many years as a generator coolant with no 
serious accidents. If helium is used as the coolant, a 
marked penalty is paid in blower power requirement; 
this calculates to be in excess of a 100° increase. Leakage 
of the hydrogen/oxygen into the atmosphere/coolant 
systems represents the greatest hazard, and it is partly for 
this reason (as well as for methane and ozone removal) 
that the platinum catalyst chamber was added to con- 
vert these gases to water and carbon dioxide. Short of 
major destruction of some part of the system, the use of 
hydrogen as a heat transfer agent should present no 
danger. Detection equipments exist which could act, in 
case of excess leakage, to dump the hydrogen and replace 
it from the helium tank. This would allow time for 
descent from orbit at any selected point along the orbit 
path. 


Vv. CONCLUSION 

As shown in Table II, the weight of the system designed 
for a 60-day manned satellite is relatively low when 
compared with other systems. Some developmental 
research is required to make it operable, but only minor 
changes should be necessary to obtain a satisfactory 
system. While designed for one man, the addition of 
automatic controls would make it suitable for an animal- 
carrying satellite. Further, it can be readily scaled for 
any crew size desired. 

For trips longer than about three months, the use of a 
calcium oxide or equivalent regenerative system for 
carbon dioxide removal is definitely indicated. Such a 
system requires considerable power and heat dissipation, 
as the temperature of most carbonate decompositions is 
about 1000°C. 

It appears reasonable to assume that from reasons of 
safety helium as a heat exchange gas would be more 
attractive for a manned vehicle than hydrogen. How- 
ever, as indicated in the Appendix, a mass and power 
penalty is involved in its use. If mass is the deciding 
factor, then hydrogen is apparently the most efficient 
heat transfer substance available. 

When the vehicle descends from orbit, the hydrogen 
coolant is dumped overboard along with the radiator 
system. During the descent, cabin cooling is main- 
tained by evaporation of the water reserve. It is assumed 
that the thermal flux into the cabin from re-entry heating 
is reasonably low. 


APPENDIX 

For long-term operation, radiation cooling is the only 
feasible method of removing heat from a body in a 
vacuum. The familiar Stefan-Boltzmann law (E = 
GeT*) allows a simple solution for a radiator whose 
emissivity is known. As deep space has a radiation 
temperature of about 5°K., it can be assumed to be a 
perfect sink for all but very low black body tempera- 
tures.* 

The temperatures assumed for the heat exchange gas 
are inlet 80°C., outlet —S5°C. As the radiator must be 
at a lower temperature than the gas to absorb heat, it was 
necessary to arrive at some temperature which would 
average the radiator black body temperature. By 
treating the surface in infinitesimals, an average tempera- 
ture of ~30°C. was arrived at. Due to unknowns, this 
value is likely off by + 10%. However, sufficient allow- 
ance has been made to cover such errors. 

An emissivity « of 0-8 was chosen as this is readily 
achievable. Using a G of 0-482 x 10-" B.T.U./ft.*/ 
sec./(°R.)*, the average radiation level was found to 
be 0-482 x 107% x 5:86 x 10” x 0-8 = 2:26 x 10° 
B.T.U./ft.2/sec. 

The B.T.U. input from all sources is assumed as 
150,000 B.T.U./24 hr. However, as most of this arrives 
during the sunlight period, a basic radiation capacity of 
2:5 B.T.U./sec. is assumed necessary. This requires a 
surface radiating area of approximately 110 ft.2 Adding 
30% to this to allow for loss of gas channels due to 
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meteor puncture or mechanical leaks gives a total surface 
area of approximately 142 ft.* 

The radiator material chosen is black surfaced alu- 
minium alloy approximately 0-01 in. thick with an 
overall mass of ~0-13 Ib./ft.2 This gives a total skin 
mass of ~37 ib. Ends, fittings, bladder valves and 
transducers will be about 50 Ib. additional mass, bringing 
the total to 87 Ib. 

To transfer the peak of 2:5 B.T.U./sec. to the radiator 
at an inlet temperature of 80°C. and an outlet of —5°C. 
requires a flow rate at constant volume at | atm. of 
0-31 Ib./min. hydrogen or 1-27 Ib./min. helium. For 
the stated conditions and 30°C. average skin tempera- 
ture, this is 59 ft. hydrogen or 127 ft.* helium. Air 
ft.3/min. x 46/in.? x 144 

33,000 ‘ 
at 0-5 lb. pressure head this is 59 x 0-5 x 0-435 x 10°? 
== 0-128 H.P., and for helium 127 x 0-5 x 0-435 x 10°? 
= 0-276 H.P. Assuming a 50% efficient impeller and a 
50%, efficient motor, a peak power requirement of 380 W. 
= 0-276 H.P. Assuming a 50% efficient impeller and a 
50% efficient motor, a peak power requirement of 
380 W. for hydrogen or 825 W. for helium is required. 
However, based on the average heat load of 1-74 B.T.U./ 
sec., the average blower power requirement is 265 W. for 
hydrogen and 575 W. for helium. To insure proper 
physical contact between the gas and radiator inner 
surface, it is necessary that turbulent flow exist. This, in 
turn, slightly raises the blower power requirement due to 
increased back pressure. As the gas must remain in 





H.P. is given by For hydrogen 


physical contact with the radiator skin for a time sufficient 
to transfer its heat, it is necessary for the radiator to 
have some minimum gas volume. For hydrogen this is 
~6 ft. and for helium ~10 ft. The volume available in 
specified flow channels is adequate for hydrogen but 
leaves no margin for helium. The channels are designed 
with a |-in. internal wall-to-wall spacing and to occupy 
80% of the radiator surface area. 

Compressor power to operate the bladder valves 
calculates out on the basis of 1 ft.* of adiabatic 10 Ib./in.? 
(gauge) gas per hour as 0-01 H.P. A constant 10 W. 
consumption is assumed. A compressor and tankage 
mass of 4 Ib. is assumed as obtainable for this system. 
8 lb. mass is assumed for the valves, tubes, and Ledex- 
type stepper which controls the gas flow to and from 
the bladder valves. The Ledex valve relay will use 
~100 W. per step. 
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CORRESPONDENCE 


(In these letters to the Editor, some passages have necessarily been omitted) 


PROFAC Power Sources: a Correction 


Sir, 

The following two corrections should be made to my 
paper’. In Section II.5 (p. 119), third paragraph, the 
fourth sentence should read: “Such a power source with 
auxiliary equipment would weigh about 42 tons with the 
present state of the art and perhaps as low as 20 tons with 
the expected development of nuclear power sources in 
10 years... .” 

The correct masses are 42 tons and 20 tons, not 1] 
tons and 2 tons as given in the paper. 

Yours very truly, 
S. T. DEMETRIADES. 
Northrop Corporation, NORAIR Division, 


1001 East Broadway, Hawthorne, 
California, U.S.A. 


17 November, 1959. 


REFERENCE 
1. S. T. Demetriades, J.B.1.S., 1959-60, 17, 114. 


Earth Satellite Problem’ 
Sir, 

It is required to place a given payload into a circulat 
orbit of a specified radius by means of a two-stage “all- 
propellent” rocket with minimum propellent consumption. 
It is assumed that the Earth is non-rotating and spheri- 
cally symmetrical ; that both stages burn instantaneously, 
one at the surface of the Earth and the second at satellite 
altitude ; and that the flight is in vacuum. The payload 
ratio (ratio of total propellent mass to payload mass) is 
expressed by: 


R=e .. i xt (1) 


where 


1 Jtl s 
F=— ,/=| V21—p,\1—p2) + pV | 
rey (1—p,(1—p.*) + peVp (2) 
and ¢ = tan~ [V(p,/p.)* — 1] 

¢ being the direction angle from horizontal of first-stage 
thrust (launch direction). 
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2 
@ = tn (32) = 1 
where 
o,= 2 _ ry fadius of 
earth = 3960 miles 
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C\—First stage exit velocity (miles / sec) 
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C2 —Second stage exit velocity (miles /sec) 


Fic. 2. Nomogram for optimum direction angle ¢ from 
horizontal of first-stage thrust (launch direction) for the case 
Pe SPr 


The two nomograms reproduced in Fig. | and 2 solve 
the above equations quickly and with a sufficient degree 
of accuracy. 

A. S. LEVENS, 
Professor of Mechanical Engineering. 


Mechanical Engineering Mechanics and Design, 
University of California, 
Berkeley 4, California, U.S.A. 


14 October, 1959. 


REFERENCE 
1. G. Leitmann, J.B.1.S., 1959-60, 17, 194. 


Manned Navigation and Guidance in the Solar 
System 


Sir, 

The interesting paper by Cashmore and Gordon! 
calls for some comment on the subject of accuracy 
required in calculations or measurements of position 
and velocity. There are several distinct problems here, 
and an astronomer may perhaps be forgiven for the 
comment that the authors seem to be unduly pessimistic 
about the ones which we do know something about, and 
surprisingly optimistic about the others! 

The essential requirement in any form of orbit deter- 
mination is a knowledge of velocity at a given place and 
time. The shape, size and orientation of the orbit 


depend solely on the velocity; unless this is known with 
some precision the orbit is badly determined and no 
predictions of the future course of the vehicle are of any 
value. There are, of course, two different methods for 
obtaining these velocities; we can either difference 
successive positions of the vehicle, or we can integrate 
successive values of the accelerations. Both schemes 
are precisely the same in theory, but can give very diff- 
erent results in practice. The differencing scheme leads 
at once to one of the standard methods used for finding 
an orbit from three observations of one body at well- 
spaced intervals of time. Unfortunately the method is 
of little value if the apparent path of the observed body 
lies on a great circle, and this is almost certain to be the 
case with any observations of one planet made from a 
vehicle moving almost in the ecliptic. 

However, a long series of observations, followed by 
a least-squares solution, would give continuously better 
results, but would still fall short of the ideal. If we 
assume with the authors that we can fix positions within 
1000 miles, then velocities determined by differencing 
will be of the order of 1000 miles/day in error. It is 
true that the accuracy will increase with the number of 
fixes, but it is still not good enough to make predictions 
for months ahead. It may also be remarked that 1000 
miles/day (roughly 0-01 mile/sec.) is hardly likely to be 
reached by any measures of stellar aberration. There 
is, perhaps, some promise in the idea of using the radio 
Doppler shift, but this is for the future, when there are 
many radio beacons, each with fully stabilized oscil- 
lators—and a stabilized receiver in the space vehicle. 

The concept of an orbit, with its angular measures of 
node, anomaly, etc., is very much out of date, and is 
quite unnecessary. It lingers on in the literature because 
it affords a mental picture of the orbit in space, but a 
much more efficient method of representing the instan- 
taneous orbit of a body in space is to quote its three 
coordinates and three velocity components at a given 
time. With these it is a simple matter to integrate the 
motion for future times, and this method is standard 
practice in dealing with comets and minor planets; the 
attractions of the Sun and planets along three axes are 
integrated to give velocities, and these integrated again 
to give positions. The process is precisely the same as 
that used for dealing with accelerometer reacings, 
although the units may be different. If the initial data 
are sufficiently accurate the integration may be pushed 
as far ahead as is necessary, and it is quite usual for a 
complete ephemeris of a minor planet for a year ahead 
to be produced on a large digital computer in about 
10 minutes. 

The routine is, in fact, so simple that a large number of 
precomputed orbits could be prepared in a short time 
and stored on tape or film. If we combine these inte- 
gration schemes with those obtained from the accelero- 
meters, we have a complete representation of the path 
of the vehicle in space. Here, however, we come to a 
major difficulty, for the two parts of the scheme are of 
very different orders of magnitude. The methods for 
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dealing with the gravitational forces are well established 
and are not as difficult as the authors suggest. The 
forces are all quite small, and the largest of them—the 
attraction of the Sun—is only about g/2000 at the dis- 
tance of the Earth. It would be quite sufficient, for an 
interplanetary journey of the type visualized in this paper, 
to have the attractions of the Sun, and Jupiter (and 
possibly Saturn) only, and to have these at an interval 
of ten days. Interpolation at this interval is simple; 
there is certainly no need for the interval of 25 seconds 
suggested by the authors, nor for the large number of 
planets mentioned. 

The interval is not determined by the convenience of 
the user, but by the behaviour of the higher differences ; 
this is not a problem in astronomy, but of normal com- 
puting methods. The tables which we use in ordinary 
computing are designed to allow interpolation with the 
minimum mental effort, but a table of sines to ten 
decimals, at an interval of 6°, may readily be interpolated 
if the user is willing to go out to eighth differences. A 
digital computer can do this with ease, but if the interval 
is too large the higher differences become difficult to 
handle, and may eventually be divergent, when inter- 
polation is impossible. 

It is for this reason that accelerometer readings are 
difficult to assimilate into the main scheme, for here the 
forces concerned are very large, and are applied for a 
very short time, so that the rate of change is enormous. 
It is not generally realised how large these thrusts can 
be. Thus, a force of 10 g applied for precisely one 
second will cause an increase in velocity of 218 mile/hr. 
(0-061 mile/sec.) and an increase in distance of over 5000 
miles by the end of a day. And this does not take into 
account the fact that the size and shape of the orbit 
have been changed, and that the thrust, if not applied 
tangentially, will also change the whole orientation of 
the orbit in space. Such large and sudden forces could 
not possibly be entered into the integration scheme of 
the gravitational field; there must be two schemes with 
widely different intervals, and the two sets of results 
could then be added together by standard methods at 
fixed intervals. It is much more important that the 
accelerometers should be able to measure the three 
components of thrust continuously (rather than in steps) 
over very short time-ranges, and that they should be able 
to do so very accurately. A constant error of only 
0-00001 g in the readings will give a velocity error of 
455 mile/day after only one day, and a positional error 
of over two million miles after 100 days. 

This suggests that no single method is sufficient. 
Perhaps a combination of the two will be useful—con- 
tinuous fixes which can be used as a check on the inte- 
gration schemes. But this is only part of the story. 
Suppose the navigator finds himself on the wrong course ; 
what does he do then? We are told that he then 
calculates the changes in velocity which he must apply, 
but how is this done? It may be possible to derive 
algebraic equations which show that it can be done, 
but this is a very different matter, as any computer 


knows, from being able to do it. Quite frankly, I can 
think of no method, other than trial and error, by which 
the required three components of velocity can be calcu- 
lated. Here is a case where a large number of pre- 
computed orbits would be of value; they take up very 
little room, and it would not be difficult to devise a simple 
differential correction which would convert the present 
orbit into an improved one which would get the vehicle 
to about the right place at about the right time. The 
navigator then has to apply a thrust in a direction at a 
computed angle to the course of the ship for a limited 
time, but as he has only his (approximate) integration 
tables to give him that course, which cannot otherwise 
be measured, the correction also will be a mere approxi- 
mation. 

Let us face up to the fact that the methods of navi- 
gating in space that have been devised up to the present 
are about equivalent to those used by the old Elizabethan 
sailors; the navigator of a space vehicle in free fall will 
need just as much enterprise and must be just as willing 
to take risks as were the old sea-dogs of Drake’s day. 
The solution to the problem is to evolve methods which 
do not demand frequent corrections of the course and 
here the ideal of constant drive is obviously the one at 
which to aim. 

Yours faithfully, 
J. G. PORTER. 
Royal Greenwich Observatory, 
Herstmonceux Castle, 
Hailsham, Sussex. 


REFERENCE 
1. D.J. Cashmore and C. N. Gordon, J.B.1.S., 1959-60, 17, 46. 


Sir, 


We have read with interest the letters from Professor 
Lawden,' Dr. Noton,? and Dr. Porter* which comment 
on our paper.* Since some of the points made by these 
correspondents tend to overlap, we propose to reply to 
all three in this letter. 

We are grateful to both Professor Lawden and Dr. 
Noton for having our attention called to the two excellent 
papers to which they refer.** There is undoubtedly a 
great deal to be said for not forcing the vehicle back on 
to the precalculated reference trajectory, especially when 
fuel is at a premium. The linear perturbation methods 
of correction set out in these two references are very 
similar, and it is interesting to note that they are exten- 
sions of a guidance system which could be used for 
ballistic missiles. The flexibility of such a system de- 
pends to some extent upon the range of position and 
velocity parameters over which the linear perturbation 
theory is accurate enough. This touches upon Professor 
Lawden’s second point concerning the right balance 
between guidance accuracy and computation complexity. 
Perhaps a useful criterion would be that the initial 
guidance accuracy should be such that linear perturba- 
tion techniques could subsequently be used without 
significant error. 








384 CORRESPONDENCE 


However, as Dr. Porter points out in the penultimate 
paragraph of his letter, the navigational process can 
only be one of successive approximation. The require- 
ment on each phase of navigation are set by the next. 
The boost phase must be good enough for the linear 
perturbation theory to be sufficiently accurate in the 
mid-course phase. Similarly, the mid-course navigation 
and guidance must be accurate enough for the terminal 
guidance to become operative early enough with injec- 
tion errors small enough for the vehicle to be able to 
correct. 

We must agree with Dr. Noton’s first point concerning 
radio guidance. In the early stages of space exploration, 
however, it is not realistic to expect radio tracking devices 
to be operative on the other planets, nor can we expect 
Earth-based radio systems to give adequate accuracy 
at very long ranges, for example take-off from a planet 
for a return journey. It seems wise therefore to make 
the navigator primarily self-contained, but able to make 
use of extra aids such as radio wherever it is available. 

Dr. Noton may be right as to the approximation for 
computing V and D. This is basically a question of 
accuracy versus complexity of equipment and this can 
only be judged by a detailed systems assessment for 
each specific project. 

The difficulty with stellar aberration is that even along 
the reference trajectory the velocity vector changes in 
magnitude and direction and therefore the distortion of 
the celestial sphere will vary with time. If the angular 
errors (perhaps as large as 1’ arc) are not acceptable 
during the periods that navigational corrections are 
being applied by the motors, a possible solution is to 
insert biases into the electrical outputs of the tracking 
telescopes. These biases might be precomputed on the 
basis of the reference trajectory, stored in the computer 
and changed as simple functions of time. 

Again, Dr. Noton may be right on the question of the 
type of computer to be used. It is dangerous to be 
dogmatic at this early stage and only considerably more 
thought and experience can settle such matters. Dr. 
Porter’s point concerning the different problems involved 
in the integration of thrust and gravitational accelera- 
tions is an example of a factor bearing on the nature of 
the computer. It would certainly appear that two sep- 
arate integration processes are required here. In fact 
it would seem that a complete astro-intertial navigator 
feeding into a separate dead-reckoning computer is 
called for. 

We were glad to see Dr. Porter dispensing with the 
classical concepts of orbits. The straightforward method 
of computing in x, y, z co-ordinates is perhaps less 
graphic but is certainly more suited to inertial navigation 
and computing. 

We agree entirely with Dr. Porter that there is a real 
problem in determining velocity. This is perhaps the 
central problem in space navigation as we can foresee it 
today. At what point do position fixes give better 
velocity information than that available from dead 
reckoning? Clearly a very careful study must be made 


of the balance of payments between increasing the 
accuracy and frequency of planetary fixes and their 
reduction on the one hand and increasing the complexity 
and accuracy of the inertial navigator (and any other 
system) used in the boost phase and the dead reckoning 
computer on the other. 
Yours truly, 
D. J. CASHMORE. 
C. N. GorDon. 


86 Marston Gardens, Luton, Bedfordshire. 
5 Swiftsgreen Road, Luton, Bedfordshire. 


20 April, 1960. 
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Meteoritic Dust and Rainfall 


Sir, 

On suggestion of Mr. Alexander Fucs, Chairman of 
the Associagdo Brasileira de Astronomia (A.B.A.) I 
had the satisfaction to give a resumé before this Society 
in Rio de Janeiro of the excellent and very interesting 
paper “Meteoritic Dust and Ground Simulation of 
Impact on Space Vehicles” by Mr. D. H. Robey.'! Dur- 
ing our discussion of this paper were pointed out some 
essential points, that may perhaps interest you, the author 
and eventually also your readers. Therefore I have 
been charged to communicate them, as follows: 

(1) As shown in the annexed Table I, extracted from 
my worldwide collection of daily rainfall figures, 
characteristic Brazilian Stations show rainfall 
peaks exactly on 17 December, 1955, correspond- 
ing to Bowen’s Rainfall Correlation with the 
Leonid meteor shower. If allowance is made for 
a tolerance of 1 to 3 days of the rainfall effect 
from the exact date, as it appears in accordance 
with the mechanism involved (and as Bowen 
admits), then also the atmospheric precipitation 
peaks immediately following the exact date in 
Washington (U.S.A.), San Francisco (U.S.A.), 
Ziirich-Kloten (Switzerland) and Sydney (Aus- 
tralia) may be attributed to the meteoritic dust 
of the Leonid meteor shower. The daily rainfall 
figures at the last December dates seem to be 
related with Bielid I meteor shower. 

(2) Mr. Robey’s new explanation of Bowen’s Rain- 
fall Correlation together with the alluded measure- 
ment of the density of micrometeorites in the 
high atmosphere by the Aerobee Rocket, confirm 
the sound fundament of Bowen’s theory, on which 
I have developed the “Universal Meteor-Astro- 
nomical Laws of Atmospherical Precipitation.”*-* 
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TaBLeE 1.—Some Examples of Daily Figures of Atmospheric Precipitation about the Date of 17 December, 1955 








December: 10 11 12 13 14 #15 6 17 18 1s 20 21 a a oo fe ee ee ee Fee 

Rio pve Janeiro D.F. (Jardin 

Botanico), Brazil at on 1 — 3438 06 - 18-8 o2 — — — 98 29 —- -— 108 14158 75 4mm. 
FORTALEZA (Ceara), Brazil —- - - —J/ 26 wit? 2o - —- — — — 20 10 10 10 24 — mm. 
Morro VELHO (Minas Sues 

Brazil .. _—- — — 229 30 236 89 406254 180 03 03 56 150 109 1:5 191 399 Si 490 351 432 mm. 
WASHINGTON, U.S. A. ne —- — - — 0-03 — - — _-_ — — - — Tet — — in. 
San Francisco (Cal.), U.S.A. .. — — T 006 002 T 099 1:20 - 075 2-99 1-59 023 O15 1-33 016 T T 036 005 in. 
ZiricH-KLOTEN (Switzerland) .. 19 10 5 1 6 8 5 5 emg = i — il 3 a 6 — 2 _ mm. 
Sypney (Australia) ee _- 19 1 0 2 _ - 105 135 197 10 — 8 — 10 39 15 137 33 Pointst 





*T = Traces. t 1 point = 0-01 in. = 0-254 mm. 


(3) Discussing already in 1956 in the World Meteoro- 
logical Organization in Geneva the extremely 
world-wide spread of meteoritic condensation 
nuclei even to antipode stations to the direction 
of influx of the meteor streams, according to 
world-wide daily rainfall statistics, I suggested 
the Earth’s gravitational attraction force to be 
responsible for this phenomenon. In the mean- 
time the recent satellite experiments have shown 
in practice that objects in satellite orbits may 
revolve several times a day round the Earth. 
This, as well as Mr. Robey’s conclusion (p. 29) 
that meteoroid dust could remain longer in 
satellite orbits, seem to confirm my above- 
mentioned suggestion. 

(4) Mr. Robey mentions (in Table II, p. 29 of his 
paper) “allergy” as a phenomenon in connection 
with freezing nuclei. We should be interested 
in any proved evidence of this kind. 


Yours sincerely, 
HuGo MELZER. 
Icarai Palace Hotel, Praia das Flechas, 
Niteroi, E. do Rio, Brazil. 
10 July, 1959. 
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Electrical Theory of Gravity 
Sir, 


I wish to thank Mr. Cashmore and Mr. Gordon for 
bringing up some fundamental points in their letter.’ 
The “‘practical’’ unit of electric charge is the coulomb, 
which belongs to the MKS system of units wherein the 
coulomb is defined, not by Coulomb’s law, but in terms 
of the ampere, the unit of electrical current used in all 
scientific and engineering measurements.* As no test 
charge is involved in the definition of the ampere, the 
existing definition of the electric charge accepted above 
would need no remodelling to take account of the 
electrical theory of gravity presented in my paper.* 

However, the coulomb may be expressed in terms of 
the statcoulomb, which is defined by Coulomb’s law as 
as that charge which repels an equal charge of the same 
sign with a force of one dyne when the charges are 


separated by one centimetre. The precise relationship? 
is: 


2-99592 x 10° statcoulombs. 


Now let us consider what would happen to this 
relationship if we modified the definition of statcoulomb 
in terms of the attraction between two equal charges of 
opposite electrical sign. Then we find that this force 
of attraction would be | part in 10* larger than the 
force of repulsion (between electrical charges of equal 
magnitude and like sign) to account for the electrical 
theory of gravity. This means that the “modified 
statcoulomb” (for unlike charges) should be | part in 
10'* smaller than the present statcoulomb defined above. 
Or to put the other way around, if we were to remodel 
the definition of statcoulomb in terms of a force of 
attraction of one dyne between unlike charges separated 
by one centimetre, then the present statcoulomb would be 
1-000 000 000 000 000 001 of this “modified statcou- 
lomb.” This figure, for all present precision experi- 
ments is the same as unity and the relationship between 
the “modified statcoulomb” and the coulomb would be 


1 coulomb = 2-99592 x 10° “modified statcoulombs” 


which is the same as that for the statcoulomb to the 
precision that we are concerned with here or even to the 
eighteenth place of decimals if ever we can achieve that 
degree of precision. Hence since there is so little differ- 
ence between the “modified statcoulomb” and the 
present statcoulomb there would not be any necessity 
to remodel the definition of the statcoulomb either. 

The unit electric field will also not be much affected 
by the electrical theory of gravity. The basic definition 
of unit electric field in the MKS system is that a force 
having a magnitude of | newton would be exerted on a 
test charge of | coulomb placed in the field at that point. 
If the source of the above electric field and the test 
charge are of opposite electrical sign then this force 
would be one part in 10° larger than that if the source 
of the field and the test charge were of like sign (to 
account for gravity). Here again the difference is so 
small as to be negligible for all electrical measurements 
that in the interest of simplicity we may maintain the 
present definition of scalar magnitude of the unit electric 
field. The reason for trying to encourage simplicity is 
that the above definition is one of those to be first 
studied by a newcomer to electrical field theory and it 
would be imposing an additional burden for a newcomer 
to have to study two definitions instead of one, especially 


1 coulomb = 
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as none of his electrical measurements would show up 
any difference between them. 

Furthermore, maintenance of the present basic 
definition will cause no difficulty later in electromagnetic 
theory because in the calculation of the electric field in 
advanced problems, it is now the trend to compute it as 
the gradient of the potential 

Paes pa bee SEW eH 
Ox éy ez 
e.g., as in Chapter 3 of Reference (2). In such a compu- 
tation the test charge does not enter in the electric field 
and so there is no conflict with the electrical theory of 
gravity. 

Finally it appears that though there is not a 1:1 
correspondence between the unit of electric charge and 
the unit electric field, nevertheless because the electrical 
theory of gravity calls for a 


1 : 1-000 000 000 000 000 000 000 000 000 000 000 01 


correspondence which is so close to the 1:1 corre- 
spondence that for practical purposes it would be 
appropriate to keep the extensive structure of electro- 
magnetic theory as it is and treat gravity as a special 
item under electromagnetic theory and/or present the 
electrical theory of gravity in the chapter on gravitation. 

There could be very useful modifications of the ‘“‘mech- 
anisms” of electrical repulsions and electrical attractions 
so as to provide powerful propellent for spaceflight. 
But these would be the subject of patent applications and 
so for the time being would be outside the realm of pure 
electromagnetic theory. 


Yours faithfully, 


Freppy Ba HLI, 
Director of Research. 


Union of Burma Applied Research Institute, 
Kanbe, Rangoon, Burma. 


14 June, 1960. 
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Sir, 


Dr. Ba Hli’s correspondence! contains gross errors in 
physics : 


(1) The diagonal force is always smaller than the 
non-diagonal force in his Fig. 1, regardless of his 
assumption of attractive force being larger than 
the repulsive for within one part in 10°*. How 
can he justify a not attractive force? 


(2) The diagonal electric force approaches zero as 
the separation of the charges of the dipole in- 
creases to infinity. So, the net electric repulsive, 
not attractive, force varies with the moment of 


the dipole. 
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Dr. Ba Hli’s net electric force is not only not attractive 
but it is merely the result of an electrostatic couple and 
bears no relation to the unique properties of the gravi- 
tational field, any more than the Lorentz force § x B 
of a rotating body does not bear any relation to the 
gravitational field (regardless of how similar they may 
be), as some have proposed otherwise. I strongly 
recommend that Dr. Ba Hli read Vaclav Hlavaty’s 
book,? because he keeps referring to the gravitational 
field as a force. 

E. OKRess, JNR. 
81-10 135th St., Kew Gardens 35, 
New York, U.S.A. 
21 May, 1960. 
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Range and Accuracy of Ballistic Missiles: A Correction 


Sir, 

I am sorry to say that a slight error has been found in 
my paper.’ In Equation (11), p. 90, the radical should 
be preceded by a plus sign, rather than the minus sign 
shown. This does not affect any other equation, but 
Equation (11) is itself a useful final result, so the correc- 
tion is worth while. 

Yours truly, 
IRVIN G. HENRY. 
394 South Chester Avenue, 
Pasadena 5, Calif., U.S.A. 


1 May, 1960. 
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NEWS AND ANNOUNCEMENTS 


B.I.S. NEWS 


Alteration of the Bye-Laws: Officers of the Society 


At its meeting on 11 June, Council resolved that the follow- 
ing Bye-Law be adopted and added to the existing Bye-Laws 
of the Society: 

41. There shall be substituted as from 1 July, 1960, the 
titles of President and Vice-President for those of Chairman 
and Vice-Chairman respectively, provided that neither title 
shall be deemed to confer any additional powers or duties on 
the Officers concerned. The President for the time being 
shall exercise the functions of the Chairman as stated in the 
Articles of Association, and the provisions of the Bye-Laws 
with regard to Chairman and Vice-Chairman shall be read as 
if they referred to President and Vice-President respectively 
after 1 July, 1960. 

This marks a return to the pre-war practice of the Society 
and also brings the structure of the Society more into line 
with that of other professional and scientific institutions. 
When the Society was incorporated under the Companies 
Act in 1945 no provision was made for the offices of President 
and Vice-President; in recent years Council has considered 
the matter on several occasions. It is now felt that the 
change is highly desirable and will be welcomed by members 
and by other organizations with which we are associated. 

Dr. L. R. Shepherd thus becomes the first post-war President 
and the third since the Society was founded. It is particu- 
larly appropriate that he should now be the first to bear the 
title after its reintroduction, for as Chairman of Council in 
six of the last seven years he has directed our activities 
throughout a period of rapid growth and progress. 

With the adoption of the new Bye-Law, Mr. K. W. Gatland 
and Dr. W. R. Maxwell become Vice-Presidents of the 


Society. 
Scottish and Yorkshire Branches 

A special meeting of the Council was held on Saturday, 
16 July to discuss matters affecting branches and informal 
groups of the Society, and was attended by representatives 
from each United Kingdom branch. 

At this meeting two of the items on the agenda were 
applications from the Scottish and Yorkshire provisional 
branches to be granted full branch status. Council decided 
unanimously to agree to both requests, and on behalf of 
the Society the President expressed congratulations to the 
officers and members of these branches and good wishes for 
the future. 

Lecture Programme, Second Half-Session 

A number of posters advertising the Society’s lectures and 
symposia to be held during the second half of the 1960 session 
are now being printed. Members and others who have 
opportunities for displaying these on the notice boards of 
their works, colleges or offices, etc., may obtain them on 
application to the Secretary, 12 Bessborough Gardens, 
London, S.W.1. 

Posters advertising meetings of Branches are also obtain- 
able from the Branch Secretaries, whose names and addresses 
appear facing the inside of the front cover of this issue. 

It is intended that in future these posters will bear the words 


Further particulars may be obtained from . . . 


followed by a space in which the person arranging the display 
can write his or her name and address, or department, etc. 


Grant from Smithsonian Astrophysical Observatory 
The Society has received a grant of £125 from the Smith- 
sonian Astrophysical Observatory towards expenses incurred 
by Moonwatch teams during the year ended 30 June, 1960. 


This is in recognition of the useful contribution to satellite 
tracking being made by these teams. 

There are three B.I.S. Moonwatch teams, in Scotland, the 
Isle of Wight, and Bristol. Their leaders are respectively: 
Dr. A. E. Roy, Department of Astronomy, The University, 
Glasgow, W.2; Mr. S. T. Smith, Pippins Cottage, Church 
Hill, Totland Bay, Isle of Wight; and Mr. J. Carter, 7 Hur- 
lingham Road, Bristol. Members in these areas who are 
interested in participating in the work of the teams should 
write to the appropriate team leader. 

A fourth Moonwatch team is in course of formation in 
the Midlands area and interested members are invited to 
write to Mr. J. D. Llewellyn, 80 West Green Drive, Stratford- 
on-Avon, Warwickshire. 

Other teams will be established in centres where there are 
sufficient members to ensure regular manning of the observa- 
tion posts. Those interested in setting up such a team should 
write to the Secretary. 


Assistant Editor of the Journal 


Mr. A. H. S. Candlin has been appointed to a vacancy on 
the Publications Committee arising from the resignation of 
Mr. D. J. Cashmore, on grounds of ill-health. Mr. Candlin 
will also act as Assistant Editor of the Journal. 


Election of Members 
The following elections were made at the Council Meeting 
on 11 June, 1960. 


New Fellows. 


Joun ACKRoyD, B.Sc., 272, Badmi Road, D d, Bristol. 

DOUGLAS ALLANSON, 64, Kin, a Cambridge. 

PETER BAUER, Leobeth, New ‘on, Cambridge. 

“ae THOMAS BRADSHAW, AME Meche” A.F.R.Ae.S., 22, Locke King Road, 

‘eybridge, Surrey. 

J. ba ‘ONEYBEAR, M.E., Laurel Hills Road, Rt. 6, Raleigh, N.C., U.S.A. 

Epwin GerarD DAVIDSON, Byways, Kingswood Road, Tylers Green, Bucking 
shire. 

ALEXANDER GEORGE EVERARD, B.Sc., 73, New Beach Road, Darling Point, 
Sydney, Australia. 

Sin ALFRED Husert Roy Feppen, M.B.E., D.Sc., Hon. F.R.Ae.S., Hon. F.1.A.S., 
Buckland Old Mill, Bwich, Breconshire, South Wales 

STANLEY CHARLES GHOSE, M.B., ae M.A., Ph.D., DSc., M.1L.E.E., 33, Heath- 
field Road, Acton, London, W 3. 

WILLIAM Rosert Hawks, A.F.R.Ae.S., 114, Gurney Court Road, St. Albans, 
Hertf 

Epwarp W. s. Hutt, Rt. a, Many! 145, 14S, Sgetnaiaté, ' Vigne. OSs. 

Norman Eric Hype, 38, Pri 

Hersert D. JENKINS, BS., MEN as Loma Vista Place, Leo Angi, 39, Cali- 
fornia, U.S.A 

Eric KORNGOLD, B.S., M.A., 9, Albion Place, Newton Centre, Mi 

Hua Lin, A.M.LMech.E., Sc.D., 17188. 55 Street, Seattle 5, Weckingos, Sits.a. 

ARTHUR FREDERICK MCLEAN, A.M.1.Mech.E., 1604 N. Chicago Street, South Bend, 
Indiana, U.S.A. 

Cuartes Rosert SINDEN MANDeRS, M.A., Ph.D., 36, Chatsworth Gardens, 
West Harrow, Middlesex. 

EpwIn JusTIN MERRICK, M.E., McKean Road, Ambler, Penn., U.S.A. 

Georrrey Kerru CHARLES PARDOE, B.Sc., D.L.C., 2, St. James Road, Harpenden, 
Hertfordshire. 

Eric Wi_uiaM RuNDLE, A.F.R.Ae.S., 152, mg Road, London, S.E.21. 

ARTHUR HoRACE SANDFORD, 6, Sandy Close, 

CurrorD ARTHUR SEYMOUR, A.M.I.Mech.E., The British Petroleum Co. Ltd., 
B.P. House, Ropemaker Street, London, E.C.2. 

James TAYLOR, M.A., Little Frimley, Portsmouth Road, Cam , Surrey. 

Jack ALFRED Thalhimer, B.S., 6870, Meadowbrook Boulevard, Minneapolis, 26, 
Minnesota, U.S.A. 


Members elected to the Fellowship. 
, Chief ace Aerospace Division, 





Bass, 
Aeronca Mfg. Corp., Baltimore 28, Mary land, U.S.A. 
THOMAS Pease Comer, B.S., c/o owes Title Co., P.O. Box 1382, Ft. Lauder- 


dale, Florida, U.S.A. 
BARTON SEBRING PULLING, B.A., Idaho State College, Pocatello, Idaho, U.S.A. 
RICHARD STANTON Jones, “‘Ruffins,”” Pallance Road, Northwood, Cowes, Isle of 


Warten EUGENE TENGELSEN, B.S.E.E., 59, Grist Mill Lane, Halesite, New York, 
U.S.A 


New Associate Fellows. 
Josern A. ALsert, B.S.E.E., Navy BOQ, White Sands Missile Range, New 
Mexico, U.S.A. 
B. AsTuRIAs, Yepocapa, Guatemala. 
ROBERT why an “Yasme," Hillfield Road, Chalfon . 
ISTMAN, B.S., 41-07, oa Surect, ‘Emhuarst 73, New 
York, N.Y., U.S.A, 
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Basi Francis CoLiins, 72, Mt. Merrion Avenue, Belfast, 6, Northern Ireland. 

Macon C. Exuts, Jnr., B.S., Artillery Road, Yorktown, Va., U.S.A 

KENNETH HUGH FOULKEs, 58, Westbourne Grove, Yeovil, Somerset. 

Rosert Lawrence Gervais, M.S., 1000, Granville Ave., West Los Angeles 49, 
California, U.S.A. 

— GeratD Howarp, B.Sc., 41, St. Albans Road, Kingston-on-Thames, 
urrey 

RICHARD STANTON Jones, M.A., “*Ruffins,”’ Pallance Road, Cowes, Isle of Wight. 

Josern J. Kramp, 5604, Wildwood Lane, Baltimore, 9, Maryland, U.S.A 

ANTHONY JAMES Leao, B.Sc., c/o D.S.1.R. Radio Research Station, University of 
Malaya, Singapore io. 

KeiTH ALAN MaRLow, 52, Foxfield, Bandley Hill, S . Hertfordshire. 

— _— PALIN, M.A., LL.B., Avoca, St. ‘Andrew’s, Guernsey, Channel 
Siands. 

BERNARD PAtEWONSKY, M.A., M.S.E., 65, Shadybrook Lane, Princeton, New 
Jersey, U.S.A. 

a JOHN LesLie PATERSON, 18, Priory Crescent, Prittlewell, Southend-on- 

a. SSEX. 
WILLIAM zee Pe.tow, M.S., Cable & Wireless Receiving Station, 84 Ponggol 
, Singapore. 

Barry ALLEN Stein, M.A., 4711, Bedford Avenue, Brooklyn 35, N.Y., U.S.A. 

Epwin Lee Wuirte, B.E., 7, Angus Trail, N.W., = 5, Georgia, U.S.A. 

RONALD WILKES, 9, Melrose Road, Weybridge, Surre: 

PAUL VAN ZyTVELD, B.S., 6067 N. Lucerne, Chicago 46, Til., U.S.A, 





Members elected to Associate Fellowship. 


MicHaeL JaMes Scott Beton, B.Sc., ‘‘Middlefield,"’ Gainsborough, Lincolnshire. 
Awe BURTON CHATFIELD, B.Sc., 5114, Newton Street, Torrance, California, 


A. 
Joun Faust, A.B., 1516-10th Street, Berkeley 10, California, U.S.A. 
Jerry Lee Harper, M.S., General Delivery, Evergreen, Colorado, U.S.A. 
Cyrit Epwarp SHEEMAN Horsrorp, B.A., Flat 4, 26, Medway Street, London, 


Dean RUTLEDGE HULING, B.A., 2420 E. 11th Avenue Denver 6, Colorado, U.S.A. 
WituiaAM ALoysius Kinney, A.B., University Club, 1135 16th Street N.W., 
Washington 7, D.C., U.S.A 


Members elected to Senior Membership. 


Georrrey BorroMtey, 28, Merefield Street, Rochdale, Lancashire. 

ARNULF WALTER EsTerer, 932, Norwood Drive, Ashtabula, Ohio, U.S.A. 

Rosert Lionet FANTHORPE, 23, Norwich Road, Dereham, Norfolk. 

Gordon Cresswe_t Geacu, 296, Teignmouth Road, Torquay, Devon. 

RAYMOND ALBERT HARRAWAY, 38, Waldeck Street, Reading, Berkshire. 

JOHN ALEXANDER MacLeop, B.A., 4211, VanHorne Avenue, Montreal 26, 
Quebec, Canada. 

Eric ConraD Miter, 2, Helston Place, Abbots Langley, Hertfordshire. 

SrePHEN Sorsk!, 1336, Dulong, Madison Heights, Mich. , US. 

RONALD HuBert SYDNEY SPACKMAN, “Log Cabin,’ * Norwood Road, Milnerton, 
Cape Town, S. Africa. 

Wicrrep ARTHUR WuiTeHouse, 2, Church Crescent, St. Albans, Hertfordshire. 

JouN NeviLtte Woop, c/o Canadian Marconi Co., P.O. Box N. 400, St. Johns, 
Newfoundland, Canada. 


New Members. 


Donacp BassaGe, 13, Russell Gardens, Ilford, Essex. 

Sipney Georce BriGuTon, 21, Grove Park, Colingdale, London, N.W.9. 

GorDON WILLIAM BRUNT, Jeddah, Hook Hill, Freshwater, Isle of Wight. 

vom < Georce Davip Carpenter, B.Sc., 6, Sherriff Court, Sherriff Road, London, 
N.W.6. 


Joun Paut Davis, B.S., 650, Branch Street, San Luis Obispo, Calif., U.S.A. 

Daniet B. De BRA, B.E., M. S., 1421, Meadow Lane, Mt. View, California, U.S.A. 

Davip TERENCE DELANEY, 25, "Broadway Road, Bristol 7 

MArtTIN RICHARD EDMUNDS, 7, Palace View, Bromley, Kent. 

Serena Mary Fenwick, 49, Basuto Road, Parsons reen, S.W 

GeraLpD Pau Fritzke, B.S., 670, Magnolia Drive, Concord, California, U.S.A. 

Sorute Emicie Harris, 3, Park Road, Dartington, Totnes, Devon. 

WILLIAM vr vn HARTMANN, 1025, Manor Road, New Kensington, Pennsyl- 
vania, 

LYNN McNrittt Hastie, B.Sc., Christ College, Park Street, Hobart, Tasmania, 
Australia. 

GeorGce ANDRE THEODORUS HEILLEGGER, c/o Radio Caribe, Nieuwestraat 228, 
Curacao, Netherlands Antilles. 

Freperic Rocer Hearucote, 79, Reddings Road, Moseley, Birmingham 13. 

JOHN JAMES HERRINGTON, 318, Hednesford Road, Norton Canes, Nr. Cannock, 
Staffordshire. 

JoseruH Micnaet Jape, M.A., 7a, Shepherds Hill, London, N.6. 

AL = Levin, M.S.Ae., 690 Victor Way, Apt. 2, Mountain View, California, 
U.S.A. 

Davip BERNARD LOVEMAN, 462, Palmer Avenue, Teaneck, New Jersey, U.S.A. 

HuGu Rosert Lyte, The White Cottage, Kingsley Green, Nr. Haslemere, Surrey. 

Gorpon Otiver Moopy, New Zealand Electricity Department, Sub Station, 
Mataroa, New Zealand. 

ANDREW FERGUSON Nimmo, 131, Wilton Street, Glasgow, N.W. 

JosePH DesMOND Parry, B.Sc., 10, Earl Street, Keighley, Yorkshire. 

Perksins LeFevre Patton, Route 1, Box 135, Clovis, New — U.S.A. 

WILLIAM JOHN Perry, 1A, Pendoylan Street, Newtown, Cardi 

Davip Eucene Prrts, 1315 N.W. 21st, Oklahoma City 6, ee U.S.A. 

JOHN EDWARD PREEDY, 94, Greenview Avenue, Beckenham, t. 

MICHAEL QUINTON, * “Rudyard,” Lakeside Road, Lymm, Cheshire. 

Joun THomMAsS Situ, 1593 Philadelphia Street, Indiana, Pennsylvania, U.S.A. 

ALAN KLEVEN STEWART, 236, Chandler, Elmhurst, Illinois, U.S.A. 

Kerru Storey, 6, Thoresby House, Cavendish Road, Idle, Bradford, Yorkshire. 

— Gronct Stratton, B.A., Air B.P. Trading Ltd., Ropemaker Street, London, 


FRANCIS JAMES Teso, M.S., Ormlie Lodge, Thurso, Caithness, Scotland. 
ALASTAIR DE TEISSIER- PREVOST, Officers’ Mess, R.A.F. Valley, Anglesey, North 
ales. 

Georrrey Ernest THompson, M.A., Harlow Cottage, Otley Road, Harrogate, 
Yorkshire. 

BENJAMIN VAN DER Merwe, 43, Hoffe Street, Kroonstad, Orange Free State, 
South Africa. 

Renato Vesco, Via Isonso 5-12, Genova, Italy. 

DANIEL VINCENDON, 29, Boulevard Sucher, Paris, 16, France. 

CHARLES BRIAN WEATHERHOGG, “Treveth,”’ 35, Pine Road, Tividale, Dudley, 
Worcestershire. 


JOINT ACTIVITIES 


Symposium on Space Navigation 


A one-day symposium on Navigation and the Early Ex- 
ploration of the Moon is being organized jointly by the 
British Interplanetary Society and the Institute of Navigation. 
The provisional programme includes the following papers 
(some of the titles are tentative): 

“The Demands of Nearly Parabolic Lunar Trajectories on 
Differential Correction and Perturlation Theory,” by 
Professor Samuel Herrick (University of California, Los 
Angeles; a/so Principal Astrodynamicist, Aeronutronic, a 
Division of Ford Motor Company). 

“Guidance of Space Vehicles by Radio Measurements and 
Command,” by A. R. Maxwell Noton (University of 
Nottingham). 

“Role of Inertial Equipment in Lunar Flights,” by D. J. 
Cashmore (English Electric Aviation Ltd.). 

“Comparison and Implementation of Rendezvous Methods,” 
by Dr. G. R. Arthur (General Electric Company, U.S.A.). 

“Fuel Requirements for Mid-Course Corrections during 
Return from Manned Lunar Flight,” by H. L. Bloom 
(General Electric Company, U.S.A.). 

“Surface Navigation on the Moon,” by P. A. E. Stewart 
(Advanced Projects Group, Hawker Siddeley Aviation Ltd.). 

“Summary Paper,” by D. O. Fraser (English Electric Aviation 
Ltd.). 

The Symposium wil] be held in the Lecture Theatre of the 
Royal Geographical ‘society, 1 Kensington Gore, London, 
S.W.7, on Friday, 18 November, 1960. Members and others 
who wish to attend the Symposium should write to the 
Secretary of the Institute of Navigation at this address. 


European Symposium on Space Technology 


The first European Symposium on Space Technology is 
being held in London early in 1961. It is being organized 
by the Society, with the co-operation of the other principal 
European astronautical societies, and will take place in the 
Council Room of the Federation of British Industries, 21 
Tothill Street, London, S.W.1, from Monday 26 June to 
Wednesday, 28 June, 1961, inclusive. 

The Symposium is intended primarily to bring together 
scientists and engineers from Western Europe to discuss the 
possibility of formulating a joint programme of space re- 
search. A detailed programme will be given in a later issue 
of the Journal, but members and others are invited to register 
now (Registration Fee £5 5s. 0d.). Applications should be 
sent to The Secretary, British Interplanetary Society, 12 
Bessborough Gardens, London, S.W.1. 


Twelfth International Astronautical Congress 


The Twelfth International Astronautical Congress will be 
held in New York early in October, 1961, probably immedi- 
ately before the Annual Meeting of the American Rocket 
Society. It is hoped that it will be possible to hold the 
Congress in the United Nations Headquarters. 

This will be the first 1.A.F. Congress to be held outside 
Europe and should be a memorable occasion: it is possible 
that the arrangements may include a visit to the Cape Canav- 
eral launching base. It is thought that over 3000 people 
will attend the meeting; to ensure that European participa- 
tion is as large as possible, our American hosts intend to 
endeavour to arrange relatively inexpensive accommodation 
for those with limited means. 

The fare from Europe will of course be the largest item in 
the cost of attending the Congress; if sufficient people decided 
to attend, it would be possible for a special plane to be 
chartered for the occasion. If this could be done well in 
advance the fare would be of the order of £90—£100 per person. 
Members and others who think there is a possibility that they 
might attend the Twelfth Congress should let the Secretary 
know as soon as possible. 








